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The VPOS: a vast polar structure of satellite galaxies, 
globular clusters and streams around the Milky Way 
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ABSTRACT 

It has been known for a long time that the satellite galaxies of the Milky Way (MW) 
show a significant amount of phase-space correlation, they are distributed in a highly 
inclined Disc of Satellites (DoS). We have extended the previous studies on the DoS 
by analysing for the first time the orientations of streams of stars and gas, and the 
distributions of globular clusters within the halo of the MW. It is shown that the spa- 
tial distribution of MW globular clusters classified as young halo clusters (YH GC) 
is very similar to the DoS, while 7 of the 14 analysed streams align with the DoS. 
The probability to find the observed clustering of streams is only 0.3 per cent when 
assuming isotropy. The MW thus is surrounded by a vast polar structure (VPOS) of 
subsystems (satellite galaxies, globular clusters and streams), spreading from Galac- 
tocentric distances as small as 10 kpc out to 250 kpc. These findings demonstrate 
that a near-isotropic infall of cosmological sub-structure components onto the MW is 
essentially ruled out because a large number of infalling objects would have had to be 
highly correlated, to a degree not natural for dark matter sub-structures. The major- 
ity of satellites, streams and YH GCs had to be formed as a correlated population. 
This is possible in tidal tails consisting of material expelled from interacting galaxies. 
We discuss the tidal scenario for the formation of the VPOS, including successes and 
possible challenges. The potential consequences of the MW satellites being tidal dwarf 
galaxies are severe. If all the satellite galaxies and YH GCs have been formed in an 
encounter between the young MW and another gas-rich galaxy about 10-11 Gyr ago, 
then the MW does not have any luminous dark-matter substructures and the missing 
satellites problem becomes a catastrophic failure of the standard cosmological model. 

Key words: galaxies: interactions - galaxies: kinematics and dynamics - Galaxy: 
halo - Galaxy: structure - globular clusters: general - Local Group. 



1 INTRODUCTION 

That the satellite galaxies of the Milky Way fall on a great 
circle has already been put forward by iLvnden-Belll (|l976l ) 
before the advent of the standard model of cosmology. In 
recent years, this substructure defined by the 11 'classical' 
satellite galaxie^] of the MW has been confirmed and named 
the Disc of Satellites (DoS) l|Kroupa. Theis fc Boilvl 120051 : 
iMetz. Kroupa fc Jerienl 120071 ). Furthermore, the 13 fainter 
satellite galaxies , mostly de tected via the Sloan Digital Sky 
Surve y (SPSS) (lYorkll2000l). also independ ently follow this 
DoS (|Metz et al.ll2009l : iKroupa et al.ll201fj| ). The satellites 
of th e Andromeda galaxy a l so show an anisot r opic d istribu- 
tion (|Koch fc Grebelll2006l : IMetz et al.|[2007l . 120091 ), see in 
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p articular the chain-like distribution of satellites in figure 1 
of iTollerud et~aT1 (|201ll ). 

For those satellite galaxies that are fast and close 
enough such that proper motions could be determined, the 
orbital poles deri ved from their spatial kinematics support 
this anisotropy (|Metz. Kroupa fc Libeskindl 120081 ). These 
suggest that the DoS is rotationally supported as the satel- 
lites orbit within this disc. 

The satellite galaxies are popularly understood to be 
luminous dark matter halos which form abundantly in stan- 
dard cold or warm dark matter (CDM or WDM, respec- 
tively) cosmology. However, already independently of the 
spatial distributio n of the satellite galaxi es there are prob- 
lems. For example, iBovill fc Ricottll j201 ll ) point out the sig- 
nificant over-abundance of predicted b right satellites com- 
ared to observations. As shown by IStrigari fc Wechslerl 
201ll) . the bright satellite population of the MW is no statis- 
tically significant outlier compared to similar host galaxies: 
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the lack of bright satellites is common. If CDM or WDM 
cosmology were right, how can massive (~ 10 10 Mq) dark 
matter halos remain dark? 

The anisotropic distribution of MW satellites has been 
identi fied to be a challeng e for standard cosmological the- 
ory bv lKroupa et al.l (|2005l ). Subsequently, confidence in the 
CDM scenario has been sought in polar structures form- 
ing in co smological simulatio ns. This has for example been 
shown by I Brook et al.l |2008h . where torques of a neighbour- 
ing galaxy change the orientation of a galactic disc. How- 
ever, the subsequent formation of a polar disc does not re- 
solve the riddle of the MW DoS because in addition to the 
polar material, the whole population of dwarf galaxies sit- 
ting within dark ma tter subhalos would still be present. As 
iKroupa et al.l {2010) note, for satellites having i ndepe ndent 
infall and accretion histories, iLibeskind et alj (|2009T ) have 
shown that "the DM hypothesis is much less likely than 0.4 
per cent to be able to account for the MW satellite system 
in MW-type DM halos". A result that has been obtained 
by only taking into account the phase-space distribution of 
the 11 classical satellite galaxies. The probability is there- 
fore further diminished by the finding of additional systems 
aligned with the DoS. Individually infalling dark matter sub- 
halos do not align in a thin disc. The satellite galaxy system 
can therefore not be explained with any significant confi- 
dence as being dark matter subhalos with independent fall- 
in histories. 

To be nevertheless explained as dark-matter objects, 
the common phase-space distribution of the MW satellites is 
often understood as evidence of their common origin. Along 
this line of reasoning, several attempts to explain the DoS 
were put forward within the CDM paradigm, most interpret- 
ing it as a cosmological accretion structure of DM dominated 
subhalos. The satellites might have formed in a common 
group of subhalos or have fallen in along the same filament. 
These scenarios are not without problems, either. One sug- 
gestion is that the satellite galaxies fell into the MW halo as 
a group (|Li fc Helmil l20Qg| ; iD'Onghia fc Lakdl20Qg| ). How- 
ever, the DoS is too thin to be compa tible with the siz e 
of known associations of dwarf galaxies l|Metz et al.1 [2009). 
In addition, t he scenario lacks in internal consistency: on 
the one hand. [Deason et al.l (|201lt ) find that, to explain the 
coherency in velocity space, the group infall has to have 
happened recently. The only example of a group infall the 
authors show in thei r models happened at a red shift of 0.6. 
On the other hand, iNichols fc Bland-Hawthorn! (|201ll l. ig- 
noring the spatial anisotropy, find that to explain the radial 
distribution of gas-deficient and gas-rich dwarfs around the 
MW, the galaxies had to have fallen in early (at a redshift 
of 3-10). 

A second attempt (|Lovell et al.ll20lH ) suggests that fila- 
mentary accretion of the satellite galaxies results in a quasi- 
planar distribution of coh erently rotat i ng sa tellites. How- 
ever, according to fig. 1 in lLovell et al.l l|201lf ) only a maxi- 
mum of about 22 (35) per cent of the subhalofl would have 
an orbital pole within 26° (37°) of the main halo spin axis, 
which at best only barely qualifies as a coherently rotating 
satellite system. More importantly, the authors themselves 
explain that the satellites would preferentially orbit within 

2 In case of isotropy this would be 10 (20) per cent. 



the plane of the MW galaxy. In high-reso lution simu l ations 
of galaxies forming in dark matter halos, iBett et all l|201Ch 
find that it is most likely for a galaxy to have its spin vec- 
tor aligned with the direction of the angular momentum of 
the parent halo. The median misalignm ent between the two 
vectors is only « 30°. lHahn et al.l (|2010T l also find this orien- 
tation to be most likely in a simulation of galaxies forming 
in a cosmic filament, but report a larger median misalign- 
ment angle of » 50°. Th ey speculate that this discrepancy 
compared to the result of lBett et alj (|2010t ) might originate 
from resolution issues or the particular feedback recipes em- 
ployed. Another possible origin are the different large-scale 
enviro nments of the simulations, the galaxies of IBett et al.1 
(2010) mig ht reside in a l ower density environment than the 
filament of lHahn et al.1 l|2010h and have had a more qui- 
escent formation history, which might better resemble the 
quiescent merger history of the MW. 

The DoS is oriented polar to the MW disc, so almost 
perpendicular to the most probable orientation. This whole 
reasoning, however, is rendered obsolete by recalling that in 
addition to the at most 22 (35) per cent allege d coherently 
orbitin g, quasi-planar distributed DM halos in lLovell et al.l 
(|201ll ). another 78 (65) per cent subhalos are present in 
the main halo outside of this distribution, which should 
lead to about four (two) times as many satellite galax- 
ies outside the quasi-planar distribution of the DoS than 
inside. These are not observed a r ound the MW. Further- 
more, lAngus. Diaferio fc Kroupal (|201ll ). using the proper 
motions of the MW satellites, show that infall is not viable 
because the orbital angular momentum can only be dissi- 
pated through dynamical friction on the MW dark matter 
halo, implying the satellite galaxies to be hosted by unphys- 
ically massive dark matter halos. Thus, however one turns 
the coin, there is no consistent solution to the DoS problem 
in contemporary popular cosmology. 

In opposition to the dark-matter-motivated accretion 
scenario, the DoS galaxies might be of tidal origin, having 
been formed in a galaxy collision l|Kroupa et al.| [2010 V This 
scenario can be considered to be much more natural. The 
existence of the bulg e of the MW and i ts formation on a 
very short timescale (|Ballero et al.|[2007l ) already indepen- 
dently indicates that the young MW experienced an interac- 
tion, and at least partial merger, with another galaxy. The 
tidal scenario would make the MW satellite galaxies mostly 
dark-matter free ancient Tidal Dwarf Galaxies (TDGs) 
formed in an interaction of the MW progenitor and an- 
other galaxy. Their positions and orbits would then consis- 
tently lie in the plane of the i nteraction and counter-orbiting 
satell ites emerge naturally |Pawlowski. Kroupa fc de Boer] 
2011). Furthermore, in tidal interactions not only TDGs are 
produced, but also large num bers of massive star clusters 
with masses of 10 s to 1O 7 M l|Bournaud. Due fc Emselleml 
2008). These would have to be distributed around the MW 
in a similar manner. 

A purely spatial alignment of MW satellites might be 
a (however unlikely) chance result of the time of observa- 
tion. Finding a similar disc-like distribution in another set 
of supposedly accreted objects, globular clusters classified 
as young halo (YH) objects, would decrease the probability 
of the DoS being a non-genuine structure even more. In the 
tidal scenario, YH objects would have been formed together 



© 2012 RAS, MNRAS 000. [111211 



The 



VPOS: a vast polar structure around the MW 3 



with, and thus would be similarly distributed as, the satellite 
galaxies in the DoS. 

An alignment of the orbital planes of the non-satellite- 
galaxy-objects with the DoS would further support its sig- 
nificance, while a contrary result would cast doubt on the 
scenario of a pure tidal origin of the DoS. Therefore, more in- 
formation on the orbits of MW satellites is needed. One way 
to add more orbit-information is to look at the disrupted 
remnants of clusters and dwarf galaxies around the MW: 
stellar and gaseous streams. Their benefit is that the spatial 
extend of each of these structures already defines an orbital 
plane for the respective stream. Thus, no additional, usu- 
ally rather uncertain proper motion data are needed. Until 
further and more precise proper-motions for MW satellite 
galaxies are available, this method thus is a promising al- 
ternative in constraining not only the spatial, but also the 
kinematic anisotropy of the MW-halo substructure. 

The distribution of globular clusters is analysed in Sec- 
tion(2] In Section[3] at first the method for finding the stream 
normal vectors is presented (Section 13. 1|) before individual 
streams are considered (Section 13. 2|) . This is followed by a 
discussion of the preferential direction found (Section I3.3[) 
and an estimation of the probability of such a distribution 
of stream normals in case of an underlying isotropic distri- 
bution (Section I3.4|l , ending in a discussion of the effect of 
the stream-distance on the analysis (Section 13. 5[) . Section [4] 
consists of our discussion of the discovered vast polar struc- 
ture (VPOS) and the proposed scenario of the MW satellites 
being tidal dwarf galaxies. Finally, our conclusions are given 
in Section [S] 



2 GLOBULAR CLUSTER DISTRIBUTIONS 

In the following, the globular cluster (GC) system of the 
MW is analysed for a preferred spatial alignment with the 
DoS. For this, the GCs are split up into three groups. The 
old halo (OH) and bulge/disc (BD) objects are thought to 
have formed with the early MW, the latter being confined to 
the MW bulge and disc. A significant fraction of the young 
halo (YH) objects is found to be younger than t he OH clus- 
ters. They are thought to be of accretion origin (|Zinnll 19931 ; 
IParmentier et al1l2000ft . This makes them interesting in the 
context of the DoS as, when they are of similar origin non- 
native to the MW, they can be expected to show an imprint 
of this origin in their spatial distribution like the MW satel- 
lite galaxies. 

The GCs are divided into these three subsystems ac- 
cording to their metallicity and horizontal-branch mor- 
phology index. For the present analysis, the classification 
into h orizontal-branch classes by iMackev k, van den B crah 
i|2005ft is used. This results in sample sizes of 30 YH, 70 OH 
and 37 BD globular cluster s . The cluster positions are taken 
from the catalog bv lHarris! (| 19961 ) (2003 update). 

To each sample of GCs a best-fitting plane is determined 
and compare d to the DoS. The sa me method as described 
in section 5 of lKroupa et alj ((2010) to determine the normal 
vector to the best-fitting disc of the satellite galaxies is now 
applied to the GCs of the MW. In short, the GC positions 
are expressed in Galactic, Cartesian coordinates. Starting 
with a projection onto the x-z plane and rotating it around 
the z-axis (pointing to the MW north pole), for each projec- 



tion a best-fitting line is determined, the line being a plane 
seen edge-on in the current projection. The fit parameters 
of the projection with the lowest RMS height of GCs from 
the line-fit are then adopted, resulting in the direction of the 
normal to the best-fitting GC disc. 

The YH clusters individually exhibit similarities to the 
star clusters of satellite galaxies. If the YH objects are of a 
similar origin as the MW satellite galaxies, then their best- 
fitting disc should consequently show similarities to the DoS, 
while these should be absent in the OH and BD populations. 
The BD clusters however should correlate with the plane of 
the MW. The inclusion of the BD clusters thus acts as a 
consistency check, their normal has to be close to the MW 
pole. The analysis of the OH cluster sample can act as a 
test whether the found substructure is real. If observational 
biases would play a role, the OH clusters should be similarly 
affected and thus the analysis should find a similar normal 
vector for them. 

This analysis does not quantify the significance of the 
prolate-ness of the GC distributions, but merely determines 
the distributions' minor axes at each projection. To test the 
robustness of the determined normal vectors, a bootstrap- 
ping analysis is performed. For each of the three globular 
cluster subsystems, the analysis is repeated 10000 times on 
a randomly re-sampled cluster population of the same sam- 
ple size, allowing for multiple occurrences of the same clus- 
ters. The distribution of the normal vectors produced by 
this method is illustrated by the la contours in Fig. [1] Fur- 
thermore, the mean directions and spherical standard devi- 
ations from these have been determined. Finally, the shape- 
parameter 7 and strength-parameter £ dFisher et al ] |l987i ; 
iMetz et alj 12007ft allow to asses whether the distribution of 
bootstrapped normal vectors is clustered or girdled (7 > 1 
and 7 < 1, respectively), and how strongly so (the larger 
the more concentrated the distribution is, with an isotropic 
distribution having C = 0). The resulting normal vectors are 
compiled in Table ffl. 

2.1 Young Halo GCs 

For the YH clusters (blue open diamond in Fig. [1]), a strik- 
ing agreement with the DoS of the MW satellite-galaxies is 
found. The 'disc of young halo globular clusters' (DoGCyh) 
normal points towards (l,b) = (144° , —4.3°), less than 13° 
away from the DoS normal. The chance that a randomly ori- 
ented vector is this close to the DoS normal is only about 2.5 
per cent. The root-mean-square (RMS) scale height of GCs 
from the best-fitting plane is 11.8 kpc, its offset from the 
Galactic centre is 2.6 kpc. The bootstrapping analysis finds 
a mean direction of (I, b) = (145.2°, —1.3°) with standard 
deviation of 21.5° which is consistent with the DoGCyh- 
direction. The shape- and strength parameters are 7 = 3.12 
and £ = 3.36, i.e. a significant clustered distribution of boot- 
strapped normal vectors. Taken together it can be stated 

jj Note that, independe ntly and in parallel, 

iKeller. Mackev &; Da Costal l|2012h have also analysed the glob- 
ular cluster system of the MW and report similar results, even 
though their GC categorization leads to slightly different samples. 
Their full YH clusters result in a normal vector of the best-fitting 
plane pointing to (I, b) = (144° ± 6° , -13° ± 5°), while fitting the 
outer 7 YH clusters only leads to (I, b) = (162° ± 5°, -8° ± 7°). 



© 2012 RAS, MNRAS OOO.rTlfgTI 



4 Pawlowski, Pflamm-Altenburg & Kroupa 




Figure 1. Directions of normal vectors and orbital poles in Galactocentric coordinates; Imw is the Galactic longitude, Z>mw the Galactic 
latitude. Viewed from the Galactic centre, the Sun is at Zmw = 180°, &mw = 0°. For a better readability only one of the two normals for 
each object is shown (those in the range 120° < Zmw < 300°). Disc of Satellites: The filled, light magent a squares show the n ormal 
vectors for the DoS fitted to the 11 classical (vertical line) and 13 faint (horizontal line) satellite galaxies from lKroupa et al.l l|2010h . The 
fit to all 24 satellite galaxies together gives the normal vector illustrated by the open, dark magenta square. Satellite orbital poles: The 
orbital poles of the eight MW satellites analysed bv lMetz et al.l l|2008h are included as light-green dots. The solid light-green lines show 
the uncertainties in the orbital pole positions. The larger, open, green circle marks the average orbital pole of the six satellite galaxies 
that lie close to the DoS, while the thin green circle gives their spherical standard deviation of 35.4°. Globular clusters: The normal 
vector of the disc fitted to the globular clusters (see Sect. [2{ classified as young halo objects, DoGCyhi is included as the open blue 
diamond, that for the old halo clusters, DoGCohi as an open yellow diamond and that of the bulge/disc clusters, DoGCbd> as an open 
cyan diamond. The dashed contours in the respective colors represent the 1-<t regions for the normal vectors derived from bootstrapping. 
They, too, are shown in only one direction, within 120° < Zmw < 300° for the YH and OH clusters and only at the southern MW 
pole for the BD clusters. The yellow filled diamond denotes the normal vector to the OH objects when excluding the outlying cluster 
NGC 2419. The filled blue diamond with a horizontal line marks the direction of the normal vector fitted to the 20 inner YH clusters 
(Galactocentric distances less than 20 kpc), while the one with a vertical line shows the normal direction fitted to the outer YH clusters 
only (Galactocentric distances of more than 25 kpc). Streams: The filled, light-red hexagons denote the directions of the stream- normal 
vectors (see Sect. [3)1 derived from the stream anchor points compiled in Table [2] The Magellanic Stream-normal is marked by slightly 
bigger, rotated hexagon. The uncertainties derived from Monte-Carlo experiments are illustrated as the light-red areas surrounded by 
dashed 1-cr contours. The open, dark-red hexagon represents the direction of the average stream-normal, its spherical standard distance 
is illustrated with the thin red circle. There is a strong clustering close to the DoS normal, indicative of the vast polar structure (VPOS) 
around the MW. The majority of satellite galaxy orbital poles, half the stream-normal vectors and the DoGCyh normal vector all lie in 
the vicinity of the DoS normal. There is very good agreement of the DoS normal and the average stream- normal directions and they both 
are close to the mean orbital pole. The DoS normal lies within the standard deviations of the streams, orbital poles and the DoGCyh 
normal vector. The DoGCqh an d DoGCbd normals are obviously detached, with the latter aligning with the MW pole as expected. 



that the YH globular cluster distribution is indistinguish- 
able from that of the MW satellite galaxies. Both show a 
preference to lie in the same plane highly inclined to the 
MW disc. 

The fitting method is dominated by the outermost GCs. 
There is a natural radial separation point in the sample, as 
the 10 outermost YH clusters have Galactocentric distances 
of more than 25 kpc, while the inner 20 are closer than 20 
kpc. To test the domination of the 10 outermost YH clusters, 
those were fitted individually. The fit gives essentially the 



same normal vector as the full sample: (I, b) = (147°, —3.3°), 
with an angular distance from the DoS normal of now only 
9.5° and an RMS scale height of 18.3 kpc. The offset of the 
best-fitting plane for this sample from the Galactic centre is 
8.0 kpc, almo s t the same as the offset of 8.2 kpc found by 
iKroupa et "all (|201Ch for the 24 MW satellites. Both offsets 
point into the same direction (towards the reported normal 
vectors, i.e. away from the Galactic centre, approximately 
in the anti-centre direction), emphasizing the common dis- 
tribution of YH GCs and MW satellite galaxies. Interest- 
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Table 1. Directions of normal- vectors 



Name 


Symbol in Fig. [T] 


Method 




7 roi 


Objects 


Reference 


DoS 2007 


square (vertical) 


positions 


157.3 


-12.7 


11 classical satellites 


Metz et al. (2007) 


DoS 2009 


not plotted 


positions 


149.6 


-5.3 


22 satellites 


Metz et al. (2009) 


DoS 2009 (w/o Her) 


not plotted 


positions 


159.7 


-6.8 


21 satellites 


Metz et al. (2009) 


DoS 


open square 


positions 


156.4 


-2.2 


24 satellites 


Kroupa et al. (20101 


Faint DoS 


square (horizontal) 


positions 


151.4 


9.1 


13 non-classical satellites 


Kroupa et al. (20101 


Mean orbital pole 


open dot 


proper motions 


177.0 


-9.4 


6 satellites 


Metz et al. (20081 


Average stream-normal 


open hexagon 


anchor points 


155 


-26 


14 streams 


this work 


DoGCyh 


blue diamond 


positions 


144 


-4.3 


30 GC classified as YH 


this work 


DoGCyh (< 20 kpc) 


blue diamond (horizontal) 


positions 


169 


2.3 


20 inner YH GCs 


this work 


DoGCyh (> 25 kpc) 


blue diamond (vertical) 


positions 


147 


-3.3 


10 outer YH GCs 


this work 


DoGCoh 


yellow diamond 


positions 


271 


-7.4 


70 GC classified as OH 


this work 


DoGCoh (w/o NGC 2419) 


yellow diamond (vertical) 


positions 


134 


39.6 


69 GC classified as OH 


this work 


DoGCbd 


cyan diamond 


positions 


175 


-85.7 


37 GC classified as BD 


this work 



Collection of the directions of the normal vectors of planar structures around the MW, given in Galactocentric, Galactic longitude I and 
latitude b. Listed are DoS- and DoGC-fit-normals, the average satellite orbital angular momentum and the average stream-normal position. 
Directions of individual stream-normals can be found in Table [2] Note the close proximity of results obtained with independent methods 
and objects: satellite galaxy positions and proper motions, stream-normals and YH globular cluster positions follow the same VPOS. 



ingly, the planes are shifted into the approximate direction 
of the Andromeda Galaxy, which lies at galactic coordinates 
(I, b) = (121.2°, -21.6°). This is 39° away from the DoS nor- 
mal and 29° away from the normal defining the best-fitting 
plane to the outer 10 YH GCs. 

A fit to the remaining inner 20 YH GCs gives a nor- 
mal vector of (I, b) = (169°, 2.3°), 13.4° away from the DoS 
normal, with a RMS scale height of 7.3 kpc. While the two 
normals to subsets of the YH objects are 22.7° apart from 
each other, they lie on opposing sides of the DoS normal 
direction and thus are both aligned very close to the DoS. 

This suggests that the DoS is a structure that has left an 
imprint at least as close as 20 kpc from the MW centre. The 
probability to find two independent normal vectors drawn 
from a uniform distribution within 13.4° of the DoS is less 
than 0.1 per cent, a strong hint that the two YH GC samples 
share a common information in their distribution. 

2.2 Old Halo GCs 

The best fitting normal vector found for the OH clus- 
ters DoGCoh (yellow open diamond in Fig. [TJ points to 
(I, b) = (271°, -7.4°) with a RMS scale height of 4.8 kpc and 
an offset from the Galactic centre of 1.8 kpc. This, while also 
highly inclined to the MW disc, is 66° away from the DoS 
normal. However, the DoGCoh is strongly dominated by the 
outermost cluster, NGC 2419, which has a Galactocentric 
distance of about 90 kpc. The second-farthest cluster, Palo- 
mar 15, is at less than 40 kpc distance. Excluding NGC 2419, 
the resulting normal vector points to (I, b) — (134°, 39.6°), 
still 47° away from the DoS, and results in a RMS scale 
height of 3.9 kpc. The OH globular cluster distribution is 
therefore not aligned with the DoS. The bootstrapping anal- 
ysis further shows that the evidence for a best-fitting disc 
for the OH GCs is less convincing. The bootstrapped nor- 
mal vectors of the full sample show a larger scatter with 
a standard deviation of 30.7° (23.2° excluding NGC 2419). 
Their mean normal vector points to (I, b) — (279.3°, —15.8°) 
((133.1°, 40.2°) excluding NGC 2419). The angular distance 



to the normal vector of the non-bootstrapped sample of 
all OH GCs is 12°. The shape- and strength-parameters 
7 = 0.96 and £ = 2.93 (7 = 1.32 and ( = 3.62) are not con- 
clusive either, being on the edge between a clustered (7 > 1) 
and a girdled (7 < 1) distribution. Interestingly, the result- 
ing normal vector of the sample including all OH GCs points 
close to the orbital pole of the Sagittarius dwarf galaxy at 
(l,b) = (273°, -13.5°). 

2.3 Bulge/Disc GCs 

The BD clusters behave as expected. The analysis finds a 
best-fitting DoGCbd normal vector (cyan open diamond in 
Fig- H]) pointing towards (/,&) = (175°, —85.7°), very close 
to the MW pole. The RMS scale height is 1.3 kpc, the offset 
from the Galactic centre is 0.2 kpc. The bootstrapping con- 
firms this direction, with a mean bootstrapped normal vec- 
tor at (I, b) — (166.3°, —87.5°) and a standard deviation of 
only 12.7°. The shape- and strength parameters of 7 = 2.95 
and £ = 4.40 indicate further that the bootstrapped nor- 
mal vectors are strongly clustered. The BD clusters agree 
well with forming a plane co-planar to the MW disc, a re- 
sult which gives confidence in the informative value of the 
analysis. 

These findings thus support the validity of the DoS as 
a major structural component of the MW. The vastly differ- 
ing plane normal directions for the OH and the YH clusters 
and the less well-defined properties of the former confidently 
demonstrate that the agreement in normal directions of the 
DoS and the DoGCyh are not merely based on some un- 
known bias in the data. 



3 STREAMS 
3.1 Method 

The assumption that stars in tidal stre ams approximately 
follow the orbit of their parent object jOdenkirchen et al.l 
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Figure 2. Sketch illustrating the determination of the stream- 
normal vectors. The stream is illustrated by the curved, dashed 
line. The vectors s and e are pointing from the position of the 
Sun towards the start- and the end-anchor-point of the stream, re- 
spectively. They are transformed into the vectors s' and e' which 
point from the Galactic centre (GC) forwards the stream's start- 
and end-point, respectively. Their cross product results in the 
stream-normal vector n. 



2003) is supported by numerical 
globular clusters l|Combes. Leon fc Mevlan 
dSph sa tellite galaxies ( Johnsto n. Hernquist fc Boltd 



simulation s for 
1999h 



both 
and 



19961: 



IZhao. Johnston. Hernquist fc Spergell 
Montuori et al.l (|2007h found that, close to a cluster 



19991 ). 



(within 7-8 tidal radii), the tidal tail mainly point towards 
the Galactic centre and thus is a poor tracer of the cluster 
path. Longer tails, however, are good tracers of the cluster 
path, which is the case for the extended streams analysed 
in this work. For very massive satellite galaxies, it has been 
shown (|Choi, Weinberg fc Kat j [2007) that the tidal tails 
can deviate strongly from the path of the satellite. This 
deviation, however, is in the radial direction as seen from 
the Galactic centre, which does not affect our analysis of 
the stream-normal vectors because the tail remains confined 
in the orbital plane of the satellite galaxy (see their fig. 7). 

To find the stream-normal vectors, two points along the 
stream are chosen. These are, depending on the available 
data in the literature, the reported start- and end-points of 
streams or pronounced over-densities along the streams. The 
stream-normal vectors are then determined to be the normal 
vectors of a plane defined by these two anchor points and 
the centre of the MW, thus assuming the streams' centre of 
orbit is the MW centre. This assumption is warranted as the 
MW centre should be the centre of mass of the Galaxy. 

With these assumptions, the position of the normal to 
the plane of the stream is estimated. The two anchor points 
on a stream are s and e. These are two vectors pointing 
from the Sun to two parts of the stream (see Fig. [2}, given 
in Galactic, Heliocentric coordinates in Tabled In the first 
step, the vectors s and e are transformed to Galactic, Carte- 
sian coordinates. They are then shifted by <1gc, the distance 
of the Sun from the Galact i c cent re, along the x-axis. For 
consistency with lMetz et al.l (j2007j) dGC = 8.5 kpc is chosen. 
This places the origin of the Cartesian coordinate system 
for the transformed vectors s' and e' on the Galactic centre. 
The stream-normal vector n to the plane defined by the two 
points on the stream and the Galactic centre becomes 



The stream-normal vector will then be expressed in 
Galactocentric coordinates of longitude I and latitude 6, 
defining its direction. 

This method is easy to implement and can be used with 
available data for streams from the literature. The major 
sources of error are the uncertain distances to the two anchor 
points and the distance to the centre of the MW. Further- 
more, in cases of extremely diffuse streams spreading over a 
large angular scale, the angular position of the anchor points 
can be badly defined. To estimate these uncertainties, we 
adopt a Monte Carlo method. We vary the anchor-point an- 
gular positions and radial distances according to a Gaussian 
distribution with FWHM of the uncertainties as compiled 
in Table [2] The resulting uncertainty in the position of the 
stream-normal vector for 10000 Monte-Carlo realisations is 
illustrated by 1-er contours in Fig. [T] 

For each stream-normal vector, the angular distance 
from the D oS normal vector p ointing to I — 156.4° and 
b = -2.2° dKroupa et al.1 12010| ) is calculated. The uncer- 
tainties in this angular distance are not necessarily symmet- 
ric. An extreme example is the stream of NGC 5053, which 
has uncertainty-contours that are almost tangential to the 
DoS normal and thus more extended to higher angular dis- 
tances from the DoS normal than to lower. Therefore, the 
uncertainties are expressed as follows: the lower (upper) un- 
certainty interval is that interval within which fall ler of the 
Monte-Carlo realisations that lead to smaller (larger) dis- 
tances from the DoS normal than the non-varied case. The 
results are listed in Table [2] 

3.2 Individual stream data 

The method described is applied to 14 streams (stellar and 
gaseous) in the MW halo for which anchor points with an- 
gular distances of more than 5° could be extracted from the 
literature. No shorter streams or mere tidal deformations 
were taken into account. The individual streams, the choice 
of anchor points and possible problems or difficulties are dis- 
cussed below. The stream data used are compiled in Table El 
together with the respective references. In total, all known_| 
(and suspected) stream-structures close to the MW are in- 
cluded in the analysis. Reasons for leaving out the Virgo 
Stellar Stream (VSS) from the present analysis are given at 
the end of this section. 

The streams can be divided into two groups, one having 
stream-normal vectors pointing close to the DoS normal (less 
than 35° away) and those further away (more than 50°). In 
the following, they are sorted by their stream-normal vec- 
tor's angular distance from the DoS normal, starting with 
the closest stream. 

GCN: The GCN high velo c ity clo ud complex shows 
a string of HI clouds which Ijinl (|201dh show to probably 
be a gaseous stream. The author fits an orbit to a number 
of GCN clouds. Jin's initial-condition point for the orbit 
estimation is taken as one anchor point for the stream. Its 
radial distance results from the orbit calculation. The second 
anchor point is the position of a compact high velocity cloud 

4 To the authors. 
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Table 2. Streams of the Milky Way 



Name 


Type 


z HC H 


6 HC [0] 


Alb[°] 


r HC [kpc] 


8[°] 




bg C [°] 


e[°] 


ref 


GCN stream 


GS 


39.3 
18 


-31.0 
47 


±1 
±3 


25 ± 5 
11 ± 1 


80 


161 


-14 


13±I 


(7) 


NGC 5053 stream 


GC 


312.7 
335.6 


77.7 
79.0 


±0 
±1 


16.4 
16.4 ± 1.64 


5 


158 


-16 


14tf 


(8) 


Styx 


DG 


311.5 
42.7 


82.8 
30.0 


±1.7 
±1.7 


38.0 ±7.6 
50 ± 10 


60 


140 


-3 


17 ±3 


(3) 


Magellanic Stream 


DG/GS 


73.4 
29.4 


-68.0 
-83.4 


2.0 
±2.0 


55.0 ± 11.0 
55.0 ± 11.0 


IS 


179 


3 


24±l° 


(1) 


Cocytos 


GC 


289.5 
41.6 


59.2 
29.6 


±1 
±1 


11 ± 2 
11 ± 2 


75 


165 


-28 


27tl 


(3) 


Pal 5 stream 


GC 


355.3 
22.3 


45.9 
36.3 


±0.5 
±0.5 


23.2 ± 2.0 
23.2 ± 2.0 


22 


163 


25 


27+_t 


(4) 


Lethe 


GC 


232.7 
41.2 


67.9 
30.5 


±0.5 
±0.5 


12.2 ± 2.4 
13.4 ±2.6 


81 


145 


-32 


32± 5 3 


(3) 


GD-1 


GC 


196.7 
113.9 


47.4 
58.1 


±0.25 
±0.25 


7.7 ±0.5 
9.1 ± 0.5 


48 


129 


-52.3 


55 ±1 


(5) 


Orphan 


DG 


249.5 
173 


50.0 
46.5 


±0.7 
±0.7 


21.4 ± 1.0 
46.8 ±4.5 


49 


203 


-47 


60±l 


(11) 


NGC 5466 stream 


GC 


42.2 
154.9 


73.6 
71.8 


±0.0 
±2 


16.3 ±0.0 
16.3 ± 3.3 


29 


276 


13 


60 ±4 


(6) 


Sagittarius 


DG 












273.8 


-13.5 


64 


(9) 


Anticentre 


DG 


224.8 
151.0 


20.5 
37.6 


±2.5 
±2.5 


8.9 ± 0.2 
8.9 ± 0.2 


65 


143 


-68 


67 ± 1 


(2) 


Acheron 


GC 


359.8 
42.7 


43.9 
30.0 


±0.5 
±0.5 


3.8 ±0.76 
3.5 ±0.70 


37 


230 


-55 


79 ±8 


(3) 


Cetus Polar Stream 


GC/DG 


144 
142 


-71 
-46 


±2 
±1 


36.1 ± 1.9 
30.1 ± 1.6 


25 


238 


6 


82 ±3 


(10) 



List of MW streams analysed in this work. Type lists the type or suspected origin of the stream, being either a gaseous stream (GS), a 
stellar stream of globular cluster (GC) or of dwarf galaxy (DG) origin. For each stream, the Heliocentric coordinates of two anchor points 
are listed, with Galactic longitude Z HC , Galactic latitude 6 HC , angular uncertainty A/fe[°] and Heliocentric distance r HC . These have been 
extracted from the papers listed under ref, as described in Section 13.21 The respective uncertainties are adopted from these sou rces or have 
been estimated. The Sagitarrius stream is an exception to this, its stream-normal is directly taken from[Maicwski ct al. (2003). The angle 
5 is the angular distance between the two anchor points. The resulting stream-normal vectors point to Galactic longitude Z^ c and latitude 
6§ C , this time in Galactocentric c oord inates. The angul ar di stance from the DoS normal vector is Q. 

Referen ces: flllBriins et al.N2005l); (2llGrillmair | d2006bll ; (3 )lGriIImair1 i2009ll; (4)lGrillmair fc Dionatos! d2006al l; (5 )lGrillmair fc Dionatosl 
ll2006bh: (6llGrillmair fc Johnsonl l |2006l l; (71 {20101 1; (81 lLauchner et al] ^Qul l; (91 iMaiewski et al.l l|2003l l: (101 INewberg et al.l l|2009l l~ 
(111 INewberg et al.l feoicli ! 



which Jin finds to be close to the orbit estimate in position 
and radial velocity. The radial distance to this second anchor 
po int i s estim ated from the radial distance to the orbit fitted 
bv ljinl |2010T ) (from her figs. 3 and 4), which does not vary 
much in this part of the suspected stream. 

The GCN high velocity cloud complex is the structure 
which has a stream-normal closest to the DoS normal of all 
streams analysed in this work, with an angular distance of 
only 13°. But the object is also the most speculative one. 
Comparing locations and velocities of MW satellites with 
the GCN, no possible progenitor for the stream has been 
found. It may be an ancient remnant of the gaseous tidal 
tail which formed the satellites and YH GCs. 

NGC 5053: Anal y sing SDSS data, 

lLauchner. Powell Jr. fc Wilhelml (2006) report the finding 
of a stellar tidal debris stream of approximately 1.7 kpc 
in projected total length, extending in one direction from 
the clu ster NGC 5053. This is 26 tim es the cluster's tidal 
radius jMackev fc van den Berghll2005h , therefore it can be 
expected that the tidal tail traces the orbit of the cluster. 
However, no continuous tail is detected, but a string of 



overdensities. This might indicate that the underlying 
stream is too faint to be detected such that only the 
overdensities of the tail are see n. Those might ar ise as 
epicyclic overdensities in the tail ( jKiipper et al. feoiO). The 
in dications of a tidal tail close to the cluster was confirmed 
bv lJordi fc Grebe! (|2010l l. but they did not reproduce the 
larger structure. One anchor point is chosen to be the 
cluster itself, wit h posit i on an d distance data adopted from 
the catalogue bv lHarrisl (1 19961 ). The second anchor point is 
the r eported end of the suspected tidal tail |Lauchner et al.l 
2006) assuming its Heliocentric distance to be 16.4 kpc, 
the same as the cluster. Uncertainties are only adopted 
for this second anchor point, ±1° in position and 10 per 
cent in distance. The latter results in a radial uncertainty 
similar to the projected length, the stream's orientation 
is thus varied by about ±45°. As the radial velocity of 
NGC 5053 in the Galac tic standard of rest is only 34 km/s 
IILaw fc Maiewskil |2010| ). it is reasonable to assume that 
the orbit is not aligned closely with the line of sight and 
that the tail therefore is seen more from its side than along 
its extend. 
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The stream-normal for the NGC 5053 stream points 
close to the DoS normal, with an angular distance of 14°. 
Thus, if this stream indeed follows the orbit, the cluster 
can not be associated with the Sag ittarius dwarf galaxy, 
as stated bv lLaw fc Maiewskil (|2010D . because the orbits of 
Sagittarius and NGC 5053 would be almost perpendicular 
to each other. 

Styx: The origin of this stellar stream is most likely a 
disrupted dwarf galaxy. The reported start and end positions 
|Grillmairl l2009) are adopted as anchor points, including the 
approximate distances. The width of the stream was used to 
estimate the angular uncertainty for the anchor points, the 
distance errors are estimated to be 20 per cent individually. 
A possible or igin is the Boot es III stellar over-density, also 
discovere d bvlGrillmairl ll2009h . Bootes III was later analysed 
in detail (|Carlin et alj|2009l ). leading to the conclusion that 
it is likely a currently disrupted dwarf galaxy in a transi- 
tional state b etween being a bound object and an unbound 
tidal stream (|Kroupalll997h . 

The Styx stellar stream has a stream-normal which, 
with 17° distance, is the third-closest stream to the DoS nor- 
mal. Bootes Ill's position is only « 8.5 kpc away from the 
DoS plane. The excellent agreement of the stream-normal 
(and thus the orbit of the satellite) with the normal o f the 
DoS supports therefore the claims bv lMetz et alj l|2008h that 
the DoS is rotationally supported. 

Magellanic Stream: The first anchor-point of this 
gaseous stream is chosen to be the position at which the 
two filaments of t he stream overlap, point B in fig. 12 of 
iBriins et all (|2005h . The other anchor point is point C in 
their fig. 13, marking an over-density in the Stream. The 
distance estimate to the Magellanic Stream of 55 kpc is 
adopted, deduced from the average of the distances of the 
Large and Small Magellanic Clouds (LMC & SMC). The 
distance-uncertainty for the Monte Carlo experiment is cho- 
sen to be 20 per cent, for the position-uncertainties 2° are 
adopted. 

The resulting stream-normal points close (within 15°) 
to the LMC orbital pole. As can be expected from the orbital 
poles of the LMC and SMC, the Magellanic Stream's stream- 
normal is close to the DoS, too. Its angular distance from 
the DoS normal is 24°. It is slightly offset from both the 
LMC and SMC orbital poles towards the north. 

Cocytos: This stellar stream has most likely a globular 
cluster origin. A start and an en d-position are given in the 
discovery paper (|Grillmairll2"009h . which are used as anchor 
points. The position uncertainties are estimated to be 0.5°, 
as the measured width of the stream is reported to be less 
than one degree and the positions are given to one degree 
accuracy. The distance-estimate for the stream is 11 kpc, 
with an uncertainty of 2 kpc which is adopted individually 
for both anchor points. 

Cocytos has a stream-normal close to the DoS normal, 
at an angular distance of 27°. The major part of the offset 
is in Galactic latitude, while the stream-normal is very close 
to the DoS in Galactic longitude. 

Palomar 5: this globular cluster shows rather extended 
stellar tidal tails. The anchor points are the start- and end- 
points of the stream as reported by iGrillmair fc Dionatosl 
(2006 al), as they cover the stream over a wider extent than 
lOdenkirchen et all (2001). The anchor point distances are 
adopted from their best-fit orbits and are the same as the 



distance to the cluster Palomar 5. They are allowed to vary 
by 2 kpc, while the positions are varied by 0.5°. 

Palomar 5 can be associated to the DoS as its stream- 
normal is only 27° away from the DoS normal. The cluster is 
at a distance of 14.8 kpc from the DoS plane and is classified 
as a YH object. 

Lethe: This stellar stream has most likely a globular 
cluster origin. A start and an end-positio n with distance - 
estimates are given in the discovery paper (Grillmair 2009), 
which are used as anchor points. The random distance un- 
certainties are estimated in that paper to be 10 per cent and 
the systematic one 11 per cent in case of the Acheron stream. 
As the same method is used for the Lethe distance deter- 
mination, a relative line-of-sight uncertainty of 20 per cent 
is adopted here. Note that this over-estimates the stream- 
normal vector uncertainty. The distances to the two anchor 
points are varied individually in the Monte-Carlo analysis, 
while the systematic part of the uncertainty would affect the 
distance to both anchor points in the same way. The posi- 
tional uncertainty is 0.5°, as the stream is much thinner but 
the positions are only reported to one degree accuracy. 

Lethe, like Cocytos, has a stream-normal close to the 
DoS normal, at a distance of 32° . Again, the major part of 
the offset is in Galactic latitude, while the normal is very 
close to the DoS in Galactic longitude. 

GD-1: As anchor points two positions on the stream are 
used for which IGrillmair fc Dionatosl (|2006bh give distance 
estimates. These are the midpoint of the north-eastern seg- 
ment of the stellar stream, which the authors used as the 
fiducial point for the orbit estimation, and the densest por- 
tion of the stream. The distance uncertainties are 0.5 kpc, 
the same as used by IGrillmair fc Dionatosl |2006bl ), and the 
position uncertainties for both anchor points are half the 
estimated FWHM of the stream. 

The GD-1 stream-normal is 55° away from the DoS nor- 
mal. Interestingly it is very close to the orbital plane of the 
Carina dwarf galaxy. It seems not to be associated with 
that galaxy, though, as the estimated orbital parameters 
l|Grillmair fc Dionato s 2006b) confine it to within 18 ±2 kpc, 
while the Carina dSph is at a distance of ~ 100 kpc (jMateol 
1998). Furthermore, no known GC can currently be associ- 
ated with the stream. 

Orphan: For the present analysis, the recently 
published new inves tigation of the Orphan Stream by 
INewberg et alj ( 201Ch i s used, instead of the discovery pa- 
pers (Grillmair 20063: iBelokurov et alj l2007h . Their table 
2 compiles distances and positions for different parts of in- 
creasing Galactic longitude on the stellar stream. The cho- 
sen anchor points are the first and the second-last entry in 
that table. The last entry was derived using data of only 
three stars and has t herefore not been considered, as even 
INewberg et alj (|201Cj) excluded it from their analysis. 

This stream-normal is 60° away from the DoS normal, 
but is close to the circle de scribing the stand ard deviation 
from the mean orbital pole l|Metz et alj|2008h . 

NGC 5466: There is a stellar tail associated wit h the 
globular cluster NGC 5466 l|Grillmair fc Johnsonll2006h . Un- 
fortunately, no distance estimate to the stream is given. Be- 
cause of this, the first anchor point is chosen to be the cluster 
itself , for which the distance is known to be 16.3 kpc (|Harrisl 
1996), and its position is not varied in the Monte Carlo ex- 
periment. The second anchor point is the end of the stream 
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as rep orted in the caption of fig. 1 in iGrillmair Sz Johnsonl 
(2006). For this second anchor point, the same distance as 
for the cluster is used, with a distance uncertainty of 20 per 
cent. 

While NGC 5466's distance from the DoS is only w 
2.9 kpc, its stream-normal is far from the DoS normal, at an 
angular distance of 60°. The cluster is classified as a young 
halo cluster (see Section [2]) and as such has most probably 
bee n accreted by the MW or been scattered out of the DoS 
fcf. lZhadfl99sl ). 

Sagittarius: The orbital plane o f the Sagittarius stellar 
stream was analysed in great detail (|Maiewski et al.ll2003l ; 
lFellhauer~e t al. 2006) and is in very good agreement with 
the orbital pole of the sate llite galaxy. We thus simply adopt 
the lMaiewski et all (|2003l ) stream-normal of I = 273.8°, b = 
-13.5°. 

It has already been pointed out that the Sagittarius 
satell ite galaxy is not orbiting within the DoS (|Metz et al.l 
2008), the stream-normal is at an angular distance of 64° 
from the DoS normal. However, it is close to the DoGCoh 
norm al vector. It might have been scattered on its present 
orbit <|Zhadll998l ). a possibility we discuss in Sect. 14.31 

Anticentr e: The start and end positions of the stel- 
lar stream from iGrillmairl (2006b) are used as anchor points 
in the analysis. The reported average Heliocentric distance 
and its uncertainty are adopted for both anchor points. The 
position uncertainties are estimated to be half the average 
stream width of 5°. Due to the width of the stream, the 
most likely origin is a dwarf galaxy. There seems to be sub- 
structure in the stream, with at least two more narrow com- 
ponents accompanying the wider central stream. It is spec- 
ulated that these sub-streams might be the remnants of a 
globular cluster population of the original dwarf galaxy. 

The Anticentre stream, while its stream-normal's 
Galactic longitude is close to the DoS's one, is 67° away 
in Galactic latitude. It is therefore not orbiting within the 
plane defined by the DoS. However, its orbital pole is close 
to the MW disc pole (about 20° away). 

Acheron: This stellar streams has most likely a globu- 
lar cluster origin. A start and an end-positio n with distance - 
estimates are given in the discovery paper (Grillmair 2009), 
which are used as anchor points. As in the cases of Cocytos 
and Lethe, the positions are uncertain by 0.5°. The distances 
for bo th stream positions are also reported in IGrillmairl 
(2009), their uncertainties are adopted to be 20 per cent, 
as discussed in our Section on the Lethe stream. 

The stream-normal of Acheron is, with almost 79°, the 
second-furthest away from the DoS-normal. 

Cetus Polar Stream: The anchor points for the stel- 
lar Cetus Polar Stream (CPS), are taken from the listed 
stream detections in table 1 (first and last entry) of the 
discovery paper by iNewberg. Yannv fc Willettl (|2009r i. The 
authors speculate that the CPS might be related to the glob- 
ular cluster NGC 5824, which is classified as an OH object. 

The stream-normal is the most distant one from the 
DoS normal, at an angular distance of 82°, which makes it 
almost perpendicular to the disc of satellites. 

Virgo Stellar Stream (VSS): This stellar stream or 
over-density is a problematic case and has not been included 
in the analysis as its direction is badly defined and its angu- 
lar extend i s rather short . The detection of the stream was 
reported bv lDuffau et al.l |2006T ). In that work, the stream's 



projectio n on the sky i s ali gned close to the celestial equator. 
However, iPrior et al.l (|2009l ) find four additional RR Lyrae 
stars to be likely members of the stream, which extends the 
VSS further south but does not lead to a conclusive position 
angle on the sky for the stream. More light might be shed on 
this issue by the orbit-estima te for a single RR Lyrae Sta r 
identified as a VSS member (|Casetti-Dinescu et al.l f2009^ . 
which, for reasons of consistency of the two-anchor-point 
method, was not investigated in the present work. 

3.3 Results 

In Fig.[T]the positions of the resulting stream- normal vectors 
for the analysed streams are plotted. The plot includes the 
direction of the normal t o the DoS der i ved fr om 24 satellite 
galaxies of the MW by iKroupa et all (|2010l ). In addition, 
the normal vectors determined individually for the 11 'clas- 
sical' and for the 13 'new', fa int MW satellite galaxies are 
included (|Kroupa et alJl201Ch (filled magenta squares with 
vertical and horizontal line, respectively). Also included are 
the normal vectors to the best fitting discs of globular clus- 
ters categorized as YH, OH and BD clusters (blue, yellow 
and cyan open diamon ds, respectively) and the orbital poles 
of the MW satellites (|Metz et al.ll2008l ) (green dots). Even 
though the normal vectors are two-directional since we do 
not know the orbital sense, only one direction is plotted to 
avoid a too crowded plot. Only normals within the Galac- 
tic longitude range of 120° < / < 300° are plotted. This 
range was chosen so that the Sagittarius stream-normal fits 
with the satellite's orbital pole while the majority of the 
stream-normals fall close to the average orbital pole of the 
MW satellite galaxies at (I, b) = (177°, -9°) (dark green 
diamond). 

A clustering of stream-normal vectors close to I — 160° 
is apparent, in the same direction to which the normal to 
the DoS points. Seven streams, 50 per cent of the streams 
in the analysis, are found with i n 32° of the normal vector 
to the DoS from IKroupa et all (|2010D . These are, with in- 
creasing angular distance, the GCN, NGC 5053, Styx, Mag- 
ellanic, Cocytos, Palomar 5 and Lethe streams. Thus, all 
these streams most likely orbit within the plane of the DoS. 
The anchor points of the streams aligned with the DoS have 
Galactocentric distances as low as 11 kpc. This indicates 
that there is a vast polar structure (VPOS) consisting of 
satellite galaxies, streams and GCs, which might have an 
impact on the substructure as close as about 10 kpc from the 
MW centre. The angular distance O of the stream-normal 
vectors does not correlate with the average Galactocentric 
distance of the stream anchor points. 

Spherically averaging the angular coordinates of the 
stream-normal directions of all streams results in an aver- 
age stream-normal pointing to I = 155° and b = —26°. The 
spherical standard deviation (|Metz et alj 120071 ) is osph = 
46°. This is in close proximity (24° away) to the normal 
vector of the DoS and is also close (about 27° away) to 
the average orbita l angular momentum of the MW satellites 
IjMetz et al.ll2008l ). For ease of comparison, the average satel- 
lite orbital pole, available DoS normal fits and the average 
stream-normal directions found in this work are compiled in 
Table [T] It is stunning to see how well they all agree with 
each other, given the different methods and objects used. 
In addition, the normal vector of the plane fitted to the YH 
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clusters, the DoGCyh, falls into the same region, while those 
of the OH and BD clusters are clearly distinct. 

Two of the four streams closest to the DoS most prob- 
ably have their origin in dwarf galaxies: Styx and GCN. 
Another one, the Magellanic Stream, is already known to 
emanate from the Magellanic Clouds which are known to 
orbit within the DoS. But there are also streams certainly 
or probably originating from dwarf galaxies with stream- 
normal directions further away from the DoS normal: the 
Orphan, Sagittarius, Anticentre and maybe the Cetus Polar 
Stream. 

Four of the seven streams close to the VPOS have their 
origin in globular clusters. Two of these streams are asso- 
ciated with known GCs of the MW, NGC 5053 and Palo- 
mar 5. Interestingly, both cluster are classified as YH ob- 
jects, whose distribution shows a best-fitting normal vector 
close to the DoS normal (Sect. [2}. This is an intriguing sign 
of consistency. This situation resembles the finding for the 
MW satellites, whose positions define the DoS. Those satel- 
lite galaxies with known proper motions sh ow orbital poles 
which also align with it ijMetz et alj|2008f ). Note that the 
outer clusters dominate the disc-fitting, while the two clus- 
ters with a stream aligning with the DoS are, with 16.9 and 
18.6 kpc, at intermediate Galactocentric radii for YH ob- 
jects. One third of the YH clusters have Galactocentric radii 
larger than 25 kpc. While the positions of the inner clusters 
are not crucial for the DoGCyh, their orbits might show 
signs of their accretion origin. This agrees with the finding 
that fitting the YH clusters with Galactocentric distances of 
less than 20 kpc individually also leads to a normal vector 
closely aligned with the DoS normal and thus is another sign 
that the VPOS reaches deep into the MW. 

However, the other stream associated with a YH clus- 
ter, NGC 5466 which lies at a slightly smaller Galactocen- 
tric distance of 16.2 kpc, has a stream-normal vector not 
aligned with the DoS or VPOS. It is reasonably close to 
the Sagittarius stream-normal and orbital pole, as well at 
the normal vector to t he OH globular clusters. However, 
lLaw fc Maiewskil (|2010l ') conclude that NGC 5466 could not 
have been associ ated with Sagittarius, contradi c ting e arlier 
findings by both Palma. Maiewski fc Johnston! (|2002l ) and 
iBellazzini. Ferraro fc Ibatal (|2003l ~ Scattering events might 
have changed any of the orbits though. 

The Cetus Polar Stream, with a stream-normal almost 
perpendicular to the DoS normal, might be associated with 
the OH cluster NGC 5824, which lLaw fc Maiewskil (|20101 ) 
exclude as a former Sagitarrius member. 



3.4 Probability of the stream distribution 

The clustering of stream-normals close to the a-pnori and 
independently defined DoS is intriguing. Could this be a 

chance result? 

A Kolmogorov-Smirnov test (K-S test, e.g. iPress et al.l 

1 19921 ) evaluates how likely it is that a sample is drawn from 
a given distribution function. In the present case, this null- 
hypothesis is that the streams are drawn from an isotropic 
distribution function. The angular distances, O, for the 
stream-normals from the DoS-normal have been calculated 
for each stream. Assuming an isotropic origin of the stream 
progenitors, the probability to find the observed distribu- 
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Figure 3. Cumulative probability of stream-normals within a 
fraction of the volume on the half-sphere measured from the 
DoS normal vector. Plotted are the observed stream-normal vec- 
tors (solid line) and the expectation for an isotropic distribution 
(dashed line). The two distributions clearly differ as there are 
7 out of 14 streams found within the inner 16 per cent of the 
available area. 



tion of these distances to be a random clustering can be 
estimated. 

The cumulative distribution of angular distances O of 
stream normals to the DoS normal are plotted in Fig. [3] 
The angles are expressed by which fraction of the area of 
a half-sphere is enclosed within a circle of the same radius 
as the angular distance O to the stream. In an isotropic 
distribution the probability for a stream-normal's direction 
is proportional to the area, so an isotropic distribution is 
represented by a straight, diagonal line with a slope of one 
in Fig. [3] The plot makes the reported clustering of stream- 
normals around the DoS normal obvious: 7 out of 14 streams 
are centred within less than l/6 th of the available area. 

For the K-S test, the maximum distance D between 
the cumulative probability of the stream-sample and of the 
isotropic case is calculated. In the K-S test D is a measure 
for the likelihood of the sample being drawn from the dis- 
tribution function. For the analysed sample of 14 streams, 
D is found to be 0.349. Using the approximate formula in 
IPress et al.l (|l992h . the K-S test can thus rule out the null- 
hypothesis with 95 per cent (about 2<r) confidence. 

However, the K-S test does not include the full knowl- 
edge about the clustering of the stream-normals close to the 
DoS. In fact, a similar D- value could be produced if half of 
the streams would lie far from the DoS normal, but at similar 
angular distances. The K-S test does not give the probability 
to find 7 out of 14 streams close to the DoS, but to find them 
at similar angular distances from the DoS normal. The test 
does not use the full information available, represented by 
the direction of the DoS normal vector. Therefore, a second, 
independent method to estimate the probability that the 
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observed distribution could be explained as a chance-result 
of a uniform distribution is asked for. 

It can be calculated how probable it is that at least k out 
of n stream-normals are found within an angular distance © 
to the DoS normal. The alignment of the stream-normal 
over-density with the DoS-normal has to be factored in. 
Without that, the probability that k uniformly distributed 
streams lie close to each other would be higher, but the ab- 
solute position of this clustering would be undefined. 

The probability p(G) of one vector randomly chosen on 
the sphere to point into the area Q degrees around the in- 
dependently defined DoS-normal is given by the area within 
the circle of radius O degrees on the sphere, ^4(0), in relation 
to the total area of the half-spher^fl, A(90°): 



p(e) 



A(90°) 



1 — cos O 
1 ' 



This has to be modified because the assumed isotropy 
of the normal distribution can not hold any more with 
stream-normals of high Galactic latitudes, as these are close 
to the MW poles. Streams with normals having high |6|- 
components will orbit within the MW disc, these are thus 
on the one hand hard to detect when far from the Sun, 
which will result in completeness issues, and on the other 
hand their origin could also be related to the MW disc dy- 
namics and its build-up. We therefore exclude these regions 
around the Galactic poles, defined by the angle /?, from the 
calculation of the probability, which reduces the effective to- 
tal area of the half-sphere. The probability for one normal 
vector then is given by 



p(e,/3) 



cos O 



cos {3 



Note that this is only valid if the two regions around the 
DoS normal and the MW pole do not overlap, so if the sum 
of 0, P and the distance from the DoS normal to the equator 
(2.2°) is at most 90°. 

The probability, given by a common Bernoulli experi- 
ment, for one specific set of k streams out of n being within 
O degrees to the DoS is 

P (e,/?) fc [i-p(e, p)\ n - k . 

There are, however, 



k) (n-k)\-k\ 

possibilities to draw k streams out of a total of n, allowing 
each stream to be drawn once at most independent of the 
order. Thus, the probability to find exactly k streams out of 
n within B degrees of the DoS is 



P (e,p) k [i- P (e,p)\ 



The probability to find at least k streams out of n within 
this area is then given by the sum over the probabilities to 
find k, k + 1, n streams within the area: 



5 Because the normal vectors defining the same plane can point 
in two opposing directions, so only one half-sphere has to be con- 
sidered. 




Figure 4. Sketch illustrating the minimum angle a a stream- 
normal vector can have from the axis connecting the Sun and the 
Galactic centre (GC in the sketch), which lies at Galactic coor- 
dinates (l,b) = (0°,0°). The Heliocentric distance of the stream 
anchor point is rf^eiio and the distance from the Sun to the Galac- 
tic centre is d,QQ. In this estimation, the angle a is minimal when 
the line connecting the Galactic centre is a tangent to the circle of 
radius dholio around the Sun. Thus, a stream with anchor points 
closer than c(h e ij can never have a stream normal vector closer 
than a to (l,b) = (0°,0°) and (180°, 0°). 



n— k 

E 

i=0 



k + i 



P (e,/3) ft+ Mi-p(e,/3)r 



An exclusion region of j3 = 20° around the MW pole 
makes the Anticentre stream-normal to lie barely within 
the allowed region but excludes streams orbiting within the 
disc of the MW galaxy. 7 stream-normals within © = 32° 
of the DoS are found, which has a probability of P(Q = 
32°, p = 20°, n = 14, it = 7) = 0.34 percent if the streams 
were on isotropically distributed orbits. Not considering 
the exclusion region gives a slightly smaller probability of 
P(G = 32°, P = 0°, n = 13, k = 6) = 0.24 per cent. 



3.5 Anchor point distance effects 

Fig. [4] illustrates that a stream defined by anchor points 
with Heliocentric distance dhciio smaller than the distance 
between the Sun and the Galactic centre c(gc can not lead 
to stream-normals closer to (l,b) = (180°, 0°) or (0°,0°) 
than 



a = arctan 



dac 



The DoS normal is close to (l,b) = (180°, 0°). Therefore, 
when streams close to the Sun can be assumed to have a 
higher chance of detection, streams not aligned with the DoS 
might be over-represented in the sample. 

The position of us as observers thus might prefer the 
finding of stream positions that are outside of the DoS. 
But still, the stream-normal overdensity is strong enough 
to be visible in the data analysed. The previous probability- 
estimates have thus to be considered as upper limits only. 
When excluding streams with dhelio that have no chance of 
giving normals close to DoS (only Acheron in this sample), 
the probability for the streams to be aligned in the way 
they are when drawn from a random distribution is further 
reduced to P(9 = 32°,/? = 20°, n = 13, fc = 7) = 0.20 per 
cent. 



© 2012 RAS, MNRAS 000, tWO 



12 Pawlowski, Pflamm-Altenburg & Kroupa 



VPOS edge-on 



200 


i 1 


1 

(J) o 




• 


• l' 


100 


.f 




U 

£ ° 


C 




' c 
1 ■ 

1 
1 
1 
1 
1 


'W 

-< 




1 
1 
1 


1 



' 100 200 



[kpc] 



VP05 face-on 



Q. 



200 


• 


o o 




o 


• 
• 


100 








O "^-jy 







^r- 




1 u 
' O 


o 




* \° 


■ 

■ 




• 


o 






[kpc] 


100 200 



Figure 5. The vast polar structure - VPOS - about the MW in 
Cartesian coordinates. The y-axis in both panels points towards 
the Galactic north pole. The 11 classical satellites are shown as 
large (yellow) dots, the 13 new satellites are represented by the 
smaller (green) dots, YH GCs are plotted as (blue) squares. The 
(red) curves connect the stream anchor points, the (light-red) 
shaded regions illustrating the planes defined by these and the 
Galactic centre. Note that the stream coordinates are magnified 
by a factor of 3 to ease the comparison. The obscuration-region of 
±10° from the MW disc is given by the horizontal grey areas. In 
the centre, the MW disc orientation (edge-on) is shown by a short 
horizontal (cyan) line. Upper panel: an edge-on view of the VPOS 
(x-axis aligned with I = 156.4°). The near- vertical solid line shows 
the best fit Disc of Satell ite galaxies - DoS - (seen edge-on) at 
the given projection from [Kroupa et al.l ll201Ct l. the dashed lines 
define the RMS-height of the DoS (28.9 kpc). Both YH GCs and 
streams align in a similar direction. Lower panel: a view rotated 
by 90°, the VPOS is seen face-on. The satellite galaxies and YH 
GCs show a wider distribution and the planes of most streams are 
also seen close to face-on. A movie rotating this view over 360° can 
be found online at [http: / /www. astro. uni-bonn.de/download/ 1 



4 DISCUSSION 

In this section we present our interpretation of the results 
and discuss possible explanations and implications. 

Using a new method to estimate the normal to the or- 
bital plane a stream forms around the Galactic centre, a pro- 
nounced overdensity of stream-normal vectors is unveiled. It 
is stunning that half of the streams are well aligned (having 
normal vectors closer than 32°, the next-closest one being 
55° away) with the previously known DoS, with the average 
orbital angular momentum of six nearby satellite galaxies, 
and with the normal vector to the best-fitting plane through 
the YH globular clusters. This is especially noteworthy, as 
the streams analysed not only include debris of tidally dis- 
rupted dwarf galaxies, but also consist of the remnants of 
globular clusters. The possibility of a similar alignment to 
occur by chance if the stream orbits are drawn from an 
isotropic distribution is estimated to be 0.3 per cent. 

The radial distances from the MW centre of the con- 
sidered objects show that the polar structure reaches deep 
into the MW, with stream anchor points having Galactocen- 
tric distances of 10-50 kpc, the YH GCs aligned close to the 
DoS having distances between 20 and 100 kpc and satellite 
galaxies found close to the DoS ranging as far out as 250 
kpc. As it consists of different types of objects, the struc- 
ture is termed a vast polar structure (VPOS) around the 
MW, to distinguish it from the disc of satellites (DoS) which 
only consists of the satellite dwarf galaxies of the MW. The 
three-dimensional distribution of the objects in the VPOS 
is illustrated in Fig. O in which the streams are magnified 
by a factor of three. 

4.1 Might the MW be an ancient polar-ring 
galaxy? 

The VPOS shows that the MW has a pronounced pla- 
nar, polar structure extending from 10 kpc to its outer- 
most region of 250 kpc. It can therefore be termed a polar- 
disc galaxy, a term which has also been suggested as a 
more correct replacement to the term pola r-ring galaxy (e.g. 
iGallagher et~al] |2002| ; iBrook et al] 12008ft . Could the MW 
have been a polar-ring/disc galaxy? 

The MW shares many prop erties with the criteria de- 
fined by IWhitmore et al.l jl990ft for the classification as a 
polar-ring-type galaxy. The major axes of the two compo- 
nents, the MW disc and the VPOS, are nearly orthogonal, 
as are their an gular momenta, shown by the satellite galaxy 
orbital poles (|Metz et al.ll2008ft and the stream normal vec- 
tors clustering close to the DoS normal. The centres of the 
MW and the VPOS are nearly aligned, illustrated for ex- 
ample by the small offset of only 8 kp c of the DoS from 
the Galactic centre (jKroupa et al.ll2010ft and 2.6 kpc for the 
best-fitting plane to the YH GCs. Finally, the VPOS con- 
sists of luminous material and is thin as the half-thickness or 
height of the DoS is 28.9 kpc, while the radius is ~ 250 kpc. 
Polar rings formed in numerically modelled galaxy mergers 
and interactions are stable for long times and can be trans- 
formed into polar discs when the ring is disturbed by other 
galaxies, such as t he later merger of the ring' s donor galaxy 
in an interaction jBournaud fc Combes! 12003ft . 

Speaking against the classification as a polar-ring/disc 
galaxy are the size of the VPOS (more extended than typ- 
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ical polar rings/discs) and the fact that the VPOS consists 
dominantly of individual stellar subsystems and not of more 
evenly spread-out material. However, until the long-term 
evolution of polar-ring galaxies is better understood to al- 
low a decisive answer, the MW might be interpreted as a 
candidate for an ancient remnant polar-ring galaxy. 

4.2 The polar structure in the scenario of satellite 
accretion 

The pronounced VPOS can not be expected in the accretion 
scenario of individually infalling satellite galaxies. The YH 
clusters are popularly thought to be accreted objects too, 
which might have been brought in along with dwarf galaxies, 
explaining the similar imprint in their spatial distribution. 
Both types of accreted objects will be tidally disrupted, giv- 
ing rise to streams in the MW halo. Naturally, there will be 
a preference to orbit in a similar direction, indicated by the 
stream-normal vectors aligning with the DoS and DoGCyh- 
The GCs, after having been stripped off their dwarf 
galaxy, would have orbits similar to that of their former 
host. Within the cosmological scenario, the distribution of 
YH GCs thus adds information about the preferred orbit of 
their hosts, even if those might not be around any more. 
Consequently, hosts not orbiting within the VPOS will pro- 
duce GCs not being distributed within this structure. If the 
VPOS is a young structure formed by group infall of DM 
subhalos, most YH GCs must have been stripped recently 
from dwarf galaxies belonging to that group. Then, how- 
ever, the YH GCs stripped longer ago from galaxies with 
unrelated orbits are missing. If the sa tellites fell i n alon g 
dark-matter filaments, as suggested bv lLovell et al.l (|20 1 lT ) . 
the majority of subhalos (65 to 78 per cent or more) still 
would not be distributed in the preferred plane, as their 
work has shown. As a consequence, a similar fraction of the 
dwarf galaxies populating these subhalos and therefore their 
GCs are predicted to be located and orbiting outside such a 
plane. 

In order to account for the observed large-scale align- 
ment of various subsystems within the standard cosmologi- 
cal model, individual dwarf galaxies with their GCs would 
have to fall into the MW dark matter halo in virtually the 
same plane highly inclined to the MW disc. However, there 
is no consistent explanation for a preference to accrete cos- 
mological objects in a plane perpendicular to the MW (see 
the discussion in Sect. [TJ. 

4.3 The MW satellites as TDGs 

In contrast to the cosmological scenario, a tidal scenario for 
the formation of the MW satellites as ti dal dwarf galax - 
ies (TDGs, for an introductory review see lBournaudll2010l ) 
naturally resolves the seemingly problematic correlatio n of 
satellite galaxies in phase-space l|Pawlowski et al1l201ll ). In 
case of such a tidal origin of the MW DoS, large num- 
bers of stellar cluster s form together with the TDGs, as 
iBournaud et al.l (120081 ) have shown in a high-resolution sim- 
ulation of a galaxy merger. This explains the similarity of 
the phase-space distribution of satellite galaxies, YH GCs 
and streams ranging from 250 kpc down to 10 kpc Galacto- 
centric distance in the VPOS. 



The polar structure surrounding the MW is ancient. 
Its age can be estimated from the age of its constituents. 
iMackev fc Gilmorel (|2004h have discussed the ages of the 
different popula t ions o f GCs. They used relative ages from 
ISalaris fc Weissl (|2002l ). which are given with respect to the 
age of the cluster M92. Assuming an age of 12.5 5 Gyr for 
the cluster M92 (like in lMackev fc Gilmorel 12004 ) . the ages 
of the 30 OH GCs with age estimates show a strong peak at 
about 12.2 Gyr. The 17 YH clusters with age estimates show 
two peaks, one close to the one of the OH population at 11.8 
Gyr, and a younger one at 9.3 Gyr. If the formation of the 
VPOS has been triggered by a major event, this therefore 
must have happened 9 to 12 Gyr ago. 

We suggest that the MW has experienced a near-polar 
collision with an approximately perpendicularly oriented 
disc galaxy. Figure [5] illustrates this geometry with images 
of observed interacting galaxies and two model snapshots. 
They show that polar interactions of perpendicularly ori- 
ented disc galaxies happen even at the current epoch and 
result in tidal debris distributed in a polar structure. 

The initial interaction could have been a fly-by of two 
galaxies or it might have ended in a galaxy-merger, de- 
stroying the infalling galaxy in the process. In a major 
merger, the two colliding galaxies will form a spheroidal ob- 
ject. As this event must have happened about 10 Gyr ago, 
the MW disc might h ave re-formed from gas accreted later 
l|Hammer et al.|[2005l ^ . with t he spheroidal compon ent being 
the bulge of the MW today. iBallero etail (120071 ) estimate 
that the MW bulge must have formed rapidly on a time 
scale of 0.1 Gyr, in favour of a merger-induced origin. How- 
ever, also in a fly-by encounter material stripped from the 
passing galaxy would be accreted onto the MW, possibly 
producing a bulge component. A bar instability would in 
any case channel gas onto radial orbits. 

If the TDGs were not formed in a merger but in a fly- 
by encounter, the passing galaxy still has to be nearby. Two 
candidates can be found in the Local Group: Andromeda 
and the LMC. 

The Andromeda galaxy M31 is a promising candidate 
because it is the biggest nearby galaxy and thus a large 
amount of matter could have been stripped from it to 
form TDGs. Furthermore, the dis tribution of dwa rf galaxies 
around M31 is anisotropic, too l|Metz et al.ll2007T ). with a 
disc of satellites fitted to the M31 dwarfs oriented such tha t 
it is seen nearly edge-on from the MW l|Metz et al.ll2009l ). 
This can also be seen from the positio ns of M31 dwar f galax - 
ies on the sky plotted in figure 1 of iTollerud et ail (|201lh . 
The M31 satellites are preferentially distributed in a struc- 
ture extending approximately from north to south in Galac- 
tic coordinates, just as the MW VPOS extends in north- 
south direction. A common direction of the satellite distri- 
butions of both galaxies is expected in a tidal scenario that 
formed both satellite populations together, as TDGs form 
in a plane defined by the orbit of the interaction. 

The scenario in which the satellite galaxies have been 
formed as TDGs in a past passage of the MW and M31 
therefore puts constraints on the orbits of those two galaxies, 
making it possible to predict the direction of p r oper m otion 
for M31. This was done by ISawa fc Fuiimotol ((2005), who 
suggest that the dwarf galaxies of the Local Group were 
formed in an off-centre collision of the MW and M31 some 
10 Gyr ago. While they do not consider TDG formation 
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Figure 6. Three examples of (near-) polar interactions observed at different stages. They illustrate qualitatively the possible interaction 
of the MW that resulted in the formation of its satellite galaxies as TDGs. Two approximately perpendicular disc galaxies (left picture, 
Arp 302) collide on a polar orbit. Material g ets tidally stripped and is accreted in a polar structure around the other galaxy (central 
picture, Arp 87, which Whitmo re et all lll990p ) classify as an object possibly related to polar-ring galaxies). In the tidal tails dwarf galaxies 
form (right picture, the Dentist's Chair galaxy), where the galaxies labelled A and B are the interaction partners and knots along the 
tidal tails labelled with lower-case letters are TDG candidates. They will orbit one of the two interacting galaxies in a plane defined by 
the interaction geometry. The interaction can end in a merger or in a fly-by of the two major galaxies. The lowe r row shows the particl e 
densities of two snapshots from a model of two galaxies interacting in a similar geometry (model 5deg200vel from IPawlowski ct al. 201 if) . 
They were generated by smoothing individual particle positions with a Gaussian and then plotting contours of constant density. Before 
the interaction (left) the target galaxy (red contours) is seen face-on along the plane of the interaction. The second snapshot (right) 
shows a view projected onto the plane of the interaction, with the target galaxy seen edge-on to the right. The orbital poles of model 
particles of the infalling galaxy (blue contours) at the end of this calculation are plotted in the lower panel of Figure [7] Image Credits: 
Arp 302: NASA/STScI/NR AO/ A. Evans et al ; Arp 87: NASA, ESA, and the Hubble Heritage Team (STScI/AURA); the 'Dentist's 
Chair' Galaxy AM 1353-272: IWeilbacher et al] {2003 }. 



but dwarf galaxy formation due to compression of halo gas 
during the closest approac h of the two galax ies, the orbital 
constraints are similar. As lZhao et al] (l2010t ) estimate from 
a timing argument, even in modified Newtonian dynamics 
(MOND) M31 and the MW might have had a fly-by en- 
counter in the past. Furthermore, as dynamical friction is 
much weaker in MOND (due to the absence of dark matter), 
galaxy-galaxy interactions ending in a merger are less fre- 
quent, while fly-bys can be considered more common MOND 
l|Tiret fe Combesfl20Qg| ). 

Another hint towards this scenario might be the fact 
that M31 shows a number of features which might be ex- 
plained by a violent interaction wi t h ano ther galaxy in the 
past. For example, lHammer et al] l|20ld ) have shown that 
a number of features - the morphologies and kinematical 
properties of the thin and thick disc, bulge and streams in 
the halo of M31 - can be explained by a single major merger 
about 9 Gyr ago. They also mention the formation of TDGs 
in their simulations, which might contribute to M31's satel- 
lite galaxy population or could have been eje cted into the 
direction of the MW (|Yang fc Hammed l2010l ). However, if 
only a few TDGs of M31 have been accreted by the MW, 
the probability to find them in a pre-existing, independently 
formed VPOS is low. Fouquet et al. (in prep.) therefore sug- 
gest that the 11 classical dwarf galaxies are TDGs that orig- 
inate from a major merger occurring in the history of M31. 
They formed in the first-passage tidal tail and were expelled 
into the direction of the MW where they now form a part of 



the VPOS. On the other hand, a fly-by scenario between the 
young MW and the young Andromeda galaxy should also 
be investigated. Such an interaction with an approximately 
equal-mass galaxy might have had similarly strong effects 
on the M31 galaxy and enrich both galaxies with satellite 
galaxies formed as TDGs. 

In the scenario first suggested by iLvnden-Belll (|l976l ) 
and further discussed in I Pawlowski et aU (|201ll ). the LMC 
progenitor galaxy might have been the origin of the MW 
satellite system. If the MW satellites have formed from ma- 
terial stripped from the LMC progenitor galaxy, the LMC 
needs to be on a bound orbit and at least on its second ap- 
proach to the MW. This was questioned by the recent proper 
motio n me asurements of the LMC by iKallivavalil et ahl 
l|2006l) and iPiatek et al] l|200Sl ). These revised proper mo- 
tions could resu l t in an unbound or wide o rbit |Besla et al] 
20071: IWu et all 120081: iRuzicka et ail l201Ch . which had let 
Yang fc Hammerl (|2010r ) to argue that the LMC might be 
a TDG expelled from a merger event in Andromeda. How- 
ever, one has to keep in mind that the local standard of rest 
(LSR) circular velocity is probably also higher than prev i- 
ou sly assumed jReid fc Brunthaleitl 120041 : iReid et al] I2009T I . 
As lRuzicka et al] (|2010l ) have shown, an increase in the LSR 
circular velocity increases the total mass of the MW, thus in- 
creasing the total potential energy of the LMC with respect 
to the MW, while at the same time reducing the galacto- 
centric velocity of the LMC and thus decreasing the kinetic 
energy. In total, a higher LSR circular velocity makes the 
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LMC more bound again. Consequently, assuming a higher 
LSR velocity, the Ma gellanic clouds were shown to likely 
be bound to the MW dShattow fc Loebll2009l ; iRuzicka et"aD 



l20ld ; iBartoskova et al.ll2012t ) 



Another argument in favour of a bound orbit is the con- 
nection of the LMC to the VPOS: The LMC lies within the 
DoS, the Magellanic stream shows a stream normal vector 
close to the DoS normal and the orbital pole of the LMC is 
close to the DoS normal, having an angular distance of only 
20°. This would have to be a coincidence if the LMC were 
unrelated to the VPOS and falling in for the first time. The 
probability for the LMC to show an orbital pole this close 
to the DoS normal and orbiting in the same direction as the 
majority of satellite galaxies with measured proper motions 
is only 3 per cent when assuming the probability distribu- 
tion of its infall direction on the sky is uniform. Thus it is 
unlikely that the LMC orbit is unrelated to the VPOS. 

It therefore seems to be possible that the LMC is on 
a wide but bound orbit and might have been the original 
collision partner of the MW which formed a population of 
TDGs. If the LMC is bound to the MW, this has another 
consequence. The orbits of MW satellites (GCs or galaxies) 
might be changed by scattering events, such that the present 
orbits (and consequently stream-normal vecto rs) mi g ht no t 
trace the initial configuration. In particular, IZhaol (1 19981 ) 
suggests Sagittarius to have had an encounter with the LMC 
and been scattered onto its present high-energy orbit, there- 
with solving the hitherto unsolved problem of its short or- 
bital period. If this is true, then Sagittarius might have been 
a member of the VPOS until about 2-3 Gyr ago. If the LMC 
were on its first infall, a collision of the LMC and Sagittarius 
is unlikely or impossible. While the indications for a bound 
orbit of the LMC can be interpreted in favour of such an en- 
counter, a decisive answer on the validity of the scattering 
scenario is currently not possible. 

The gaseous stream GCN, having the closest stream 
normal to the DoS, might be an ancient remnant of the 
initial tidal tail from which the satellite galaxies and globu- 
lar clusters formed. This interpretation is supported by its 
stream normal direction, which points close to the LMC or- 
bital pole and is situated in the middle between the average 
orbital pole of the satellite galaxies, the average stream nor- 
mal, the normal to the YH GCs and the DoS normal. It 
might thus indicate best t he original dir e ction of the col- 
liding galaxy. Interestingly, Wink el_ et alj (12011 ') point out 
that a model calculation of IDiaz fc Bekkil 1 201lf ) suggests 



that GCN might be associated with the leading arm of the 
Magellanic system. 



4.4 Questioning the TDG-scenario 

If the MW satellites are TDGs, then the following questions 
arise. For some, definite answers can not be given as no 
detailed studies aiming at this possible formation scenario 
for the MW satellites have been performed yet. 



4-4- 1 Do TDGs align and orbit in a disc? 

In contrast to the cosmological attempts to explain the 
VPOS, it can be considered as observationally confirmed 
that TDGs show planar distributions. As they form in tidal 





Figure 7. T he directions of th e orbital poles of the MW satellite 
galaxies from lMetz et al. (2008) and the stream normal directions 
(this work) compared to the distrib utions of orbital poles i n tidal 
debris of galaxy-interaction models ((Pawlowski ct al. 2011), plot- 
ted as contours. Note that the (-axis is shifted by 90° compared to 
Fig. [l] The upper panel shows the result of a merger calculation, 
the lower panel those of a fly-by. The orbital poles of all parti- 
cles in the models that have distances of at least 20 kpc in the 
final snapshot are included. The coordinate system of the models 
is defined in analogy to the Galactic coordinates, with the spin 
direction of the target galaxy pointing to b = —90° in the merger 
and to b = 90° in the fly-by model. The initial orbital pole of 
the infalling galaxy in both models points to (I, b) = (180°, 0°). 
The contours show the densest regions including 95, 90, 80, 70, ... 
per cent of all particles. The agreement of the observed satellite 
orbital poles and stream normal directions with the contours of 
prograde material (left, centred on I = 180°) is striking, but also 
retrograde material close to the orbital pole of Sculptor is found. 
Even the diagonal elongation of the orbital poles and stream nor- 
mals (from the upper left to the lower right) is reproduced. 



tails, observed TDGs show a strong spatial alignment with 
each other. Examples are the supermassive star clusters 
form ing in the tidal tail of the Tadpole galaxy jjarrett et alj 
2006) and the numerous TDG candidates that have been 
observed aligned alo ng the tidal tails of i nteracting galax- 
ies l ike NGC 5291 ilMalphrus et alj Il997h and AM 1353- 
272 (|Weilbacher et al.ll2002l). This spatial alignment can be 
conserved over time. iGalianni et alj d201Ch discuss that the 
satellite galaxies of NGC 1097 , embedded in a tidal stream, 
might be old TDGs. Recently. |Puc et al.l (|201ll ) found three 
old TDG candidates, aligned along a tidal tail around the 
gala xy NGC 5557. 

iPawlowski et al.l (|201ll ) have modelled collisions of 
galaxies in which a target disc galaxy, representing the early 
MW, interacts with a perpendicularly oriented disc galaxy 
falling in on a polar orbit (left plot in Figure [6} . The in- 
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Figure 8. The luminosities of MW satellites co-orbiting in the 
VPOS. The y-axis shows the V-band l uminosities Ly of the MW 
satellites in units of Lq as compiled in lKroupa et al. I fcOlCft . The 
x-axis plots the cosine of the angular distance cf> °f the respec- 
tive satellites' orbital pole from the average direction of the six 
satellites' orbital poles. For comparison, the histograms show the 
distribution of orbital poles in tidal debris in two models. They 
result from the same data as Fig. [7] assuming each particle to 
have 8,000 Lq, corresp onding to a mass-to-lig ht ratio of 2 for 
the model unit scaling in lPawIowski et"af I feOllI) . Here, cj> is mea- 
sured from the median direction of the particles' orbital poles. The 
grey error bars show the uncertainties for each bin from Poisson 
statistics. The satellites with orbital poles further away from the 
mean direction (see Table [TJ are less luminous. In the TDG sce- 
nario this might indicate that less material for their formation 
was available on these orbits, which seems to be qualitatively 
confirmed by the histograms. The amount of tidal material with 
orbital poles within a certain distance-bin follows a similar trend 
as the MW satellite luminosities. No such trend is expected if the 
satellite galaxies would have fallen onto the MW on individual 
orbits. 



teractions lead to polar discs of tidal debris, particles repre- 
senting the debris orbit within these disc and even counter- 
orbiting material is a natural occurrence, possibly explain- 
ing the counter-orbit of the Sculptor dwarf galaxy. For two 
models, a merger and a fly-by (named model 7.5degl 00vel 
and 5deg200vel respectively in iPawlowski et aL l [201ir >. the 
distributions of orbital poles of particles at the end of the 
runs is plotted in Figure [7] 

A striking agreement of the observed satellite orbital 
poles and stream normal directions with the contours of 
prograde material is apparent, even the diagonal elongation 
(from the upper left to the lower right) is reproduced. At 
lower contours, further away from the median orbital pole 
of the particles, the density of tidal material is lower. This 
might indicate that TDGs formed on such orbits might pref- 
erentially be less-luminous as less material for their forma- 
tion is available. 

The observed orbital poles of the satellite galaxies seem 
to be in agreement with this expectation, as Figure [8] illus- 
trates. The satellites with orbital poles furthest from their 
average direc tion, Carina and Draco, are the least luminous 
|Mateolll998l l. followed by the brighter Ursa Minor and For- 
nax, which are closer in, and then the the most luminous 
SMC and LMC in the centre. 

The phase-space distribution of TDGs, as expected 
from tidal debris, is therefore in good agreement with the 



observed distribution and motion of MW satellite galaxies 
in a VPOS. 



4-4-3 Do enough TDGs form in interactions? 

In many interacting galaxies young or currently forming 
star clusters and TDG candidates have b een found, for 
exam ple in the "Mice" galaxies NGC 4676 (|de Griis et al.l 
l2003h : the "Ant ennae" galaxies NGC 4038/39, NGC 3256 
and NGC 3 921 jKnierman et al l [2 0031: Arp 105, Arp 242, 



NGC 72 52 jBournaud et al.ll2004l ); t he "Tadpole" galaxy 
Arp 188 Jjarrett et al.ll2006h: Arp 305 l|Hancock et al.ll20odh 
and NGC 5557 l|Duc et aljboilf ). Numbers of TDG candi- 
dates in the order of 10 have been observed in the tidal tails 
of interacting galaxies, as in the "Den tist's Chair" galaxy 
AM 1353-272 Iweilbacher et all liool ) and in NGC 5291 
jMalphrus et al.lll997l ~ These are relatively massive object 
(typically > 10 8 Mq), equivalents of the dSph and ultra- 
faint dwarf galaxies (typically < 10 7 Mq) of the MW would 
have too low a surface brightness to be observed in distant 
galaxies. 

In a hydr odynamical high - resolu tion modelling of a 
galaxy merger, iBournaud et al.l (|2008l ) find more than 100 
new-formed objects consisting of young stars. These TDGs 
and smaller systems have masses in the range 10 5 to 10 8 Mq. 
The colliding galaxies of IBournaud et alJ (2008) are rather 
gas-poor, with a gas fraction of only 17 per cent they rep- 
resent a galaxy me rger at z ^ 1. This is important as the 
numerical models of Wetzstci n et all l|2007h have shown that 
the gas plays a major role in the formation of TDGs, such 
that an encounter in the early universe of galaxies consist- 
ing of more gas can be expected to form more or more- 
massive TDGs. The rate of TDG formation could thus be 
higher at higher redshift (|Kroupa et al.ll2010h . which seems 
to be supported by si mulations of gas-rich galaxy mergers 
IBournaud et alj|201lh . 

It it therefore plausible that enough TDGs and GCs 
form in galaxy-interactions to populate the observed VPOS 
around the MW. 



4-4-3 Can the tidal scenario explain the different 

morphologies of the LMC, the dSphs and the GCs? 

If the LMC pro genitor galaxy was the origin of the MW 



It the LMt pro genitor g alax y was the origin ot the MW 
satellite system i|Lvnden-Bell|[l976l ; [pawlowski et al.ll201ll ). 



this would explain the different morphologies of the Mag- 
ellanic Clouds compared to the dSph satellites. The LMC 
would then be the remnant of the initially infalling galaxy, 
while the majority of MW satellites would have been 
formed from the tidal debris. Even if the LMC had been 
born as a TDG, it is e xpected from numerical simulations 
ijBournaud et al.l 120081 ) that different types of stellar sys- 
tems are formed: the most massiv e, rotating TDGs a re in- 
distinguishable from dlrr g alaxies dHunter et al. 20001 ) . and 
systems below 10 8 Mq that IBournaud et al. I l|200Sft describe 
as super star clusters. These can be interpreted as com - 
pact sphero i dal ga laxies and GCs l|Bournaud et aLl feoOS). 
iBriins et al.l (|201lf ) have shown that the observed star clus- 
ter complexes forming in tidal tails of interacting galaxies 
can evolve into extended clusters and ultra-compact dwarf 
galaxies, too. 
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All this shows that a large variety of objects can form 
from tidal debris, such that the different morphologies of 
satellites observed around the MW do not rule out the tidal 
scenario of formation of the VPOS. 



4-4-4 Are the varying star formation histories of the MW 
satellites consistent with a tidal origin? 

The star formation (SF) histories of satellite galaxies born 
as TDGs can be complex and diverse. TDGs can con- 
tain a certain amount of stars from the original galaxy 
(e.g. IWetzstein et al.l 120071 ). and in addition those stars 
born when the TD G formed from the tidal-tail material. 
iRecchi et af] (|2007D have shown that ongoing SF in TDGs 
can be expected and gas ejected from the TDG might be re- 
accreted, leading to another SF period, a process that might 
be recurring. Additional gas from the tidal tail might be ac- 
creted for some time, too. Furthermore, different orbits of 
the TDGs might influence the SF rate, for example close en- 
counters with a larger galaxy can remove the gas and shut 
off SF. Therefore, complex SF histories, differing for indi- 
vidual objects born in the same interaction, are expected in 
TDGs, but no detailed studies on this are available yet. 

With the present limited knowledge about TDG forma- 
tion and evolution, the SF histories of the MW satellites do 
not provide definite hints on their origin. 



4-4-5 Can the TDGs be sufficiently long-lived? 

iBournaud fc Du e (2006) state that only 25 per cent of the 
TDGs formed in their simulations are long-lived and survive 
for more than 2 Gyr, but one needs to keep in mind that they 
only resolve TDGs with masses of at least 10 8 Mq. TDGs 
dropping below this limit are counted as non-surviving, 
while in nature less massive TDGs are found and the MW 
satellites typically are less massive, too. 

T here are ob s ervati onal hints that TDGs are stable. Re- 
cently iDuc et ail (|201ll ) found three old TDG candidates, 
aligned along a tidal tail around the galaxy NGC 5557. Their 
age estimate of these dwarf galaxies (at least 2 Gyr, but 
maybe up to 5 Gyr) demonstrates that TDGs are long-lived 
obj ects which contr i bute to the population of satellite galax- 
ies. iGalianni et~ai] (|20ld ) discuss that the satellite galaxies 
of NGC 1097, embedded in a tidal stream, might be old 
TDGs. The stability of TDGs has also been inv estigated nu- 
merica lly. Using chemodynamical calculations, IRecchi et al.l 
(2007) show that dark-matter- free dwarf galaxies easily sur- 
vive their early star formation phase, which lasts for at least 
400 Myr (no longer evolution has been modelled) . The long- 
term stabili ty of dar k -matt er-fre e satellites has been demon- 
strated by iKroupal (ll997T ) and iKlessen fc Kroupal (l998). 
Due to their lack of a massive dark matter halo, TDGs suffer 
less dynamical friction which should also help them survive 
when orbiting a more massive host. 

TDGs can therefore safely be taken to be long-lived 
objects. 



4-4-6 How can the large mass-to-light ratios of the 
satellites be explained? 

TDGs can not contain large amounts of non-baryonic 
dark matter and should therefore have lower mass-to- 
light ratios t han galaxies embedded i n massive dark - 
mat ter halos (jBarnes fc Hernquistl 1 1992ft . IKroupal (|l997f l 
and IKlessen fc Kroupa " (fT998l ) illustrated that high mass to 
light ratios of (M/L) bs > 100 can be observed in dark- 
matter free satellite galaxies, even though the intrinsic ratio 
(M/L) t ruo = 3 is low. The reasons are that the galaxies have 
non-isotropic velocity dispersions, are non-spherical and are 
not in dynamical equilibrium. This might explain the large 
observed mass-to-light ratios of some MW satellites, but it 
is improbable that all of them are out of equilibrium. 

Understanding the VPOS in terms of it being a major 
tidally formed structure is the currently only known viable 
physical process to naturally account for its phase-space cor- 
relation. This then however challenges the validity of the 
standard cosmological model (see also |Kroupa et al.l 120101: 



standard cosmological model (sec also Kroupa et al. 2U1U; 
Peebles fc Nusser 2010; Famaev fc McGaughll201ll ; IKroupal 



2012ft such that using high M/L ratios as evidence for the 



satellite galaxies being dark matter dominated would not 
be valid. For logical consistency one would need to consider 
non-Newtonian dynamics. 

It is interesting to note in this regard that even the 
observed rotation curves of young TDGs show signs of 
mi ssing mass and are su rprisingly flat, as was reported 
by IBournaud et all (120071 ). Using their data. iGentile et all 
|2007h have been able to show that the rotation curves 
are consistent with Modified Newtonian Gravity (MOND). 
In the MOND-framework, the more massive, "classi- 
cal" dSph satellites natura lly obey the Tully-Fisher rela- 
tion (|Tullv fc Fished Il977f ). as it is a basic c o nsequence 
of a MONDian gravity law llMilgroml 1 19831 ; IMcGaugh! 
120121 ; iFamaev fc McGaugbl [20111). Ind e ed, i t has a l ready 
been demonstrated dBrada fc Milgroml |2000l ; lAngusI [2008: 
iMcGaugh fc WolJ |2010| ) that the MIL ratios of the MW 
satellites would be naturally obtained in MOND, with some- 
what elevate d values for three sa tellites being the result of 
tidal effects dKroupa et al.l 120101 ; compare with th e New - 
tonian models by IKroupal ll997K IMcGaugh fc WolJ (|20ld ) 
have provided intriguing evidence that the unacceptably 
high mass-to-light ratios of the faint Local Group dwarf 
galaxies may, in the MONDian framework, result from the 
dwarfs being out of equilibrium. The preferred alignment of 
the major axes of ultra-faint dwarfs to the MW centre migh t 
hint to a tidal imprint on those satellites |Sand et al.ll201ll ) . 

Within the standard model of cosmology, current 
knowledge of TDGs does not support the observed M/L 
ratios of the MW satellites, but the issue deserves closer in- 
spection in the future. High observed M/L ratios therefore 
do not necessarily exclude a tidal origin, but can be consid- 
ered the most severe challenge of the tidal scenario if the 
dynamics are restricted to the purely Newtonian case. 



4-4-7 I s it likely that the VPOS has a polar orientation? 

Assessing the overall likelihood of obtaining a VPOS around 
a MW-like galaxy demands knowledge of numerous con- 
ditions, like those necessary for the formation of TDGs 
(a gas-rich, probably major, merger) and the probability- 
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distribution of different interaction-geometries. Many of 
these parameters in turn depend on the underlying dynam- 
ics, structure formation and environment. A decisive answer 
on the overall likelyhood can therefore not be given here, 
but some parameters are easier to evaluate than others. 

First of all, it is to be expected that about 9-12 Gyr 
ago virtually all galaxies were gas rich. If we assume that 
a galaxy-encounter leads to the formation of a structure of 
tidal debris, it then comes down to a largely geometrical 
argument whether an encounter may produce a polar struc- 
ture. For simplicity, it can be assumed that the direction of 
an infalling galaxy is drawn from an isotropic distribution 
and that the orientation of tidal debris is mostly defined by 
the orbital plane of the two interacting galaxies. In that case, 
every spacial orientation of a disc of tidally formed satellites 
is equally likely. The probability that the angle between the 
satellites disc and an arbitrarily oriented Galaxy disc is at 
least x degree is then given by cos \- Thus, the chance to be 
within 10° to a polar orientation is cos 80° = 17.4 per cent, 
or more than one in six. In half of the cases, the inclina- 
tion will be 60° or more. Thus, if a structure is formed in 
a galaxy-galaxy encounter, it is likely to be close to a polar 
orientation. 

In this regard i t migh t be interesting to note that 
iBournaud fc Combesl (I2003T ) estimate the probability for a 
galaxy to become a polar ring galaxy due to a merger to 
be at most 5 per cent. However, they also state that the 
probability to acquire a polar ring by tidal accretion in a 
fly-by encounter is significantly larger. This results in their 
prediction that many polar ring galaxies have not been de- 
tected yet as their rings are too faint. One has to exercise 
care when comparing their probabilities of polar rings to 
the VPOS. The latter structure is much larger, but also 
contains less material than a typical polar ring. One might 
expect that a larger structure is more stable as the mat- 
ter takes longer to complete one orbit and is more distant 
from the host galaxy's centre, resulting in a more symmetric 
potential, but currently this is speculation. 

While the tidal scenario can successfully explain the 
phase-space correlation of the dwarf galaxies, it is currently 
not possible to give exact numbers stating which alterna- 
tive for the formation of the MW satellite system is how 
probable. We do know however that phase-space correlated 
satellite systems are not uncommon (see Sect. 14.4. 1 p . It is 
clear that early structure formation and galaxy-galaxy en- 
counters need to be studied in order to assess the overall 
likelihood of obtaining a VPOS. But as noted in Sec. 14.4.51 
if the VPOS can most naturally be understood as a tidal 
structure, then such work needs to be performed in a non- 
Newtonian framework for logical consistency. The reason is 
that if the satellites of the MW would be ancient TDGs, 
the MW would have no DM-dominated satellites, in con- 
tradiction with the predictions of the standard cosmological 
model. 

4-4-8 Would TDGs follow the mass-size, mass-metallicity 
and other relations? 

Currently, this question can not be answered satisfactory. 
No high-resolution studies of the formation and long-time 
evolution of large numbers of TDGs have been performed. 
Trends, therefore, can not be investigated in numerical mod- 



els. Observations do not help either, as the number of known 
evolved TDGs (in contrast to actively forming ones) is still 
very low, because evolved TDGs seem to look like ordinary 
dwarf galaxies. 

If MOND (or another form of modified gravity) is valid, 
its effect on TDGs can be imagined by them being sur- 
rounded by " phantom dark matter" . This possibly results in 
similar solutions like CDM models. In addition, a stronger 
gravity would stabilize the TDGs, having the consequence 
that lifetimes and survival-rates of Newtonian dark-matter 
free models are only lower limits. 



5 CONCLUSIONS 

A vast polar structure (VPOS), consisting of satellite dwarf 
galaxies, globular clusters classified as young halo (YH) ob- 
jects, and streams of stars and gas is discovered around the 
MW. 

Applying the disc-fitting algorithm of iKroupa et al.l 

(2010) to three samples of globular clusters around the MW 
reveals a tight alignment of the normal vector to the disc 
of YH globular clusters (DoGCyh) with that of the disc of 
satellites (DoS). From this it can be concluded that the YH 
globular clusters, thought to be of extra-galactic origin, are 
distributed in a similar plane like the MW satellite galax- 
ies. These clusters and the satellite galaxies are therefore 
distributed in a common VPOS. The plane fitted to them 
is dominated by the outermost third of YH globular clus- 
ters with Galactocentric distances of more than 25 kpc, but 
fitting the inner 20 clusters independently also gives a dis- 
tribution closely aligned with the DoS. The probability to 
find the normal vectors of both the outer and the inner YH 
cluster distributions this close to the DoS normal is less than 
0.1 per cent if assuming that the normal vectors are drawn 
from a uniform distribution on the sphere. The OH clusters, 
thought to have formed with the early MW, result in a nor- 
mal vector not aligned with the DoS. Furthermore, it is much 
less certain whether their distribution can be fitted with a 
plane. As can be expected, GCs classified as Bulge/Disc ob- 
jects show a normal vector aligning with the MW disc poles, 
as they are distributed mostly in the plane of the MW. 

Using the two-anchor-point method, the normal vectors 
to streams around the MW centre have been constructed. 
These are an approximate description of a stream's orbital 
plane. It is found that half of the streams analysed (7 out 
of 14) have stream-normal vectors pointing close to the DoS 
normal vector (^ 32° away). The other half have much larger 
(^ 55°) angular distances from the DoS normal. This hints 
at a favoured orbital direction of tidal debris, in agreement 
with them preferentially orbiting within the VPOS. If the 
stream-normal vectors were drawn from an isotropic distri- 
bution, such a tight alignment with the DoS normal would 
occur in at most 0.3 per cent of the cases, ruling out this 
hypothesis with 3a confidence. Two out of three streams as- 
sociated with YH GCs are aligned with the VPOS, showing 
consistency with the GC disc fitting results, especially as 
these two clusters are relatively close to the MW and thus 
do not dominate the found DoGCyh normal vector. 

The possibility that the clustering of stream-normals is 
a random structure independent of the VPOS but by chance 
distributed in a similar direction is thus very unlikely. The 
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same holds true for the distribution of YH GCs in a similar 
plane like the DoS. In fact, for both the YH GCs and the 
streams to align with the satellite galaxy DoS has a likeli- 
hood of at most 0.1 per cent x 0.3 per cent = 3 x 10 -4 per 
cent. Therefore, both findings are significant additions to the 
DoS, not only in number but also in the type of objects that 
can be associated with it. Streams originating from dwarf 
galaxies and globular clusters alike are found close the po- 
lar structure, which turns out to be made up not only of 
satellite galaxies, but also includes globular clusters. This 
shows for the first time that not only dw arf galaxies are 
part of and orbit within (|Metz et al.l 12008? ) this structure, 
but the same is true for globular clusters. Furthermore, the 
structure is shown to range deeply into the MW, spreading 
from the farthest MW satellite galaxies at 250 kpc down to 
as close as 10 kpc from the Galactic centre for the closest 
known stream anchor points. Thus, there seems to be a co- 
herent, vast polar structure (VPOS) around the MW which 
can only have arisen in a major event including the whole 
Galaxy. 

There are s everal scenarios trying to explain the VPOS. 
A tidal origin jPawlowski et all l201ll l. in which the MW 
satellites have formed in an early, major galaxy-interaction 
from tidal debris, currently seems to be the most promis- 
ing possibility. It naturally accounts for the distribution 
of objects and their (both co- and counter-rotating) orbits 
in the plane of the interaction, which would have formed 
dwarf galaxies and star clusters together. In this scenario, 
the gaseous stream GCN may be an ancient remnant of 
the tidal tail from which satellite gal axies and star clusters 
formed. From the ages of the YH GCs jSalaris fe Weissll2002i ; 
iMackev & Gilmore!l20o4l we deduce this event to have hap- 
pened 9 to 12 Gyr ago. If the MW satellites are TDGs, we 
can now not only reconstr uct the very first format ion events 
of the MW using OH GCs jMarks fe Kroupall2010h , but also 
the one major encounter that formed the MW bulge, and 
the satellite galaxies and YH GCs distributed and orbiting 
in the VPOS. 

The alternative attempts, trying to form a major disc- 
structure around the MW from cosmological accretion only, 
suffer from insurmountable problems. Neither group-infall 
nor filamentary accretion have been shown to be able to 
account for the very pronounced, polar anisotropy in the 
spatial distribution of MW satellites in a DoS. In addition, 
they usually lead to rather diverse orbits which are not as 
strongly aligned with each other as is in fact observed. The 
increased sample-size of objects associated with the VPOS 
puts even stricter demands on their statistics, such that a 
pure chance alignment of objects resulting in the observed 
distribution can be considered ruled out. Th e implications 
of this for fundamental physics is discussed in lKroupa et al.l 
|2010l ). 
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